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ABSTRACT 


The  computerized  Quick-Reacting  General  War  Gaming  System  (QUICK)  will 
accept  input  data,  automatically  generate  global  strategic  nuclear  war 
plans,  provide  statistical  output  sunmaries  and  produce  input  tapes  to 
simulator  subsystems  external  to  QUICK.  QUICK  has  been  programed  in 
FORTRAN  for  use  on  the  CCTC  HIS  6000  computer  system. 

The  QUICK  Program  Maintenance  Manual  consists  of  four  volumes:  Volume 
I,  Data  Management  Subsystem;  Volume  II,  Weapon/Target  Identification 
Subsystem;  Volume  III,  Weapon  Allocation  Subsystem;  Volume  IV,  Sortie 
Generation  Subsystem.  The  Program  Maintenance  Manual  complements  the 
other  QUICK  Computer  System  Manuals  to  facilitate  maintenance  of  the 
war  gaming  system.^JtThis  volume,  Volume  III,  provides  the  programner/ 
analyst  with  a  technical  description  of  the  purpose,  functions,  general 
procedures,  and  programming  techniques  applicable  to  the  programs  and 
subroutines  of  the  Weapon  Allocation  subsystem.  Companion  documents  are 

A; 

a.  USERS  MANUAL 

Computer  System  Manual  UM  9-77,  Volume  I 
Computer  System  Manual  UM  9-77,  Volume  II 
Computer  System  Manual  UM  9-77,  Volume  III 
Computer  System  Manual  UM  9-77,  Volume  IV 

Provides  detailed  instructions  for  applications  of  the  system. 

b.  TECHNICAL  MEMORANDUM 
Technical  Memorandum  TM  153-77 

Provides  a  nontechnical  description  of  the  system  for  senior 
management  personnel. 


SECTION  1.  GENERAL 


1.1  Purpose 

This  volume  of  the  QUICK  Program  Maintenance  Manual  describes  the  mod¬ 
ules  which  are  part  of  the  QUICK  Weapon/Allocation  subsystem,  detailing 
the  modules,  subroutines,  and  functions  which  it  comprises.  The  infor¬ 
mation  contained  herein  is  presented  on  a  module-by-module  basis.  The 
module-by-module  discussions  are  structured  so  that  a  maintenance  pro- 
granmer  can  understand  the  program  functions  and  programming  techniques. 
The  computer  subjects  are  structured  to  inform  the  maintenance  program¬ 
mer  of  overall  system  programing  techniques  and  conventions. 

1.2  General  Description 

The  Weapon  Allocation  subsystem  operates  using  the  integrated  data  base 
as  defined  by  all  modules  up  to  PLANSET  of  the  Weapon/Target  Identifica¬ 
tion  subsystem  and  produces  a  plan  using  the  weapon  resources  specified 
to  maximize  the  expected  target  value  destroyed.  The  subsystem  consists 
of  modules  PREPALOC,  ALOC,  EVALALOC,  and  ALOCOUT,  as  shown  in  figure  1. 
Figure  2  shows  the  relationship  of  the  Weapon  Allocation  subsystem  to 
other  QUICK  subsystems  in  terms  of  procedural  and  information  flow. 

The  modules  and  supporting  subroutines  of  this  subsystem  are  used  to 
define  Information  for  use  in  later  processes  and  allocate  given  weap¬ 
ons  to  targets  to  optimize  expected  value  destroyed.  The  Integrated 
data  base  is  updated  as  each  module  is  exercised  in  sequence.  The  final 
output  provides  proper  inputs  necessary  to  execute  the  Sortie  Generation 
subsystem. 

The  first  module,  PREPALOC,  precomputes  much  of  the  information  required 
by  later  processors.  In  addition.  It  provides  capabilities  for  planning 
factor  modification  and  fixed  weapon  assignment  specification. 

The  next  module,  ALOC,  performs  the  allocation  of  weapons  to  targets. 
Using  a  generalized  Lagrange  multiplier  method,  an  optimal  allocation 
is  generated  subject  to  several  forms  of  user-input  allocation  con¬ 
straints.  These  constraints  include  specification  of  minimum  and  maxi¬ 
mum  desired  damage  levels,  restriction  of  weapons  to  specified  subsets 
of  the  target  system,  and  specification  of  weapons  allocated  to  speci¬ 
fic  targets  by  the  user.  Within  these  constraints,  the  program  gener¬ 
ates  the  allocation  which  maximizes  the  expected  value  destroyed  in  the 
target  system.  Module  ALOC  is  also  referred  to  as  the  allocator. 

The  main  function  of  module  EVALALOC  is  to  provide  a  summary  of  the 
allocation  produced  in  module  ALOC  and  to  calculate  an  expected-value 
estimate  of  its  results.  In  addition,  the  module  has  the  capability 
of  evaluating  the  effect  upon  the  results  of  variations  in  input  values 
for  weapon  and  target  parameters.  Module  EVALALOC  may  be  run  either 
before  module  ALOCOUT  or  after  module  PLANOUT. 
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Figure  1.  Major  Subsystems  of  the  QUICK  System 


2.5.3  Subroutine  SETRD.  This  subroutine  is  called  in  any  execution 
which  contains  a  SETTING  clause.  It  performs  changes  to  general  gaming 
parameters  and  weapon  group  height  of  bursts.  Changes  to  target  param¬ 
eters  are  stored  on  an  external  data  file  (25)  for  subroutine  CHGBAS . 

2.5.4  Subroutine  CHGBAS .  This  subroutine  is  called  in  any  execution 
where  subroutine  SETRD  has  stored  one  or  more  change  requests  on  file 
25.  These  requests  are  carried  out  for  the  target  base  according  to 
a  preset  priority  scheme. 

2.5C5  Subroutine  FIXWEP .  This  subroutine  is  called  in  any  execution 
with  one  or  more  FIX  clauses.  It  creates  fixed  assignment  records 
(ASSIGN)  for  use  by  the  allocation  module. 

2.5.6  Subroutine  PRNPRP ■  This  subroutine  performs  all  standard  and 
optional  prints.  It  is  called  at  the  end  of  every  run  and  produces 
those  optional  prints  requested  and  any  standard  prints  called  for  by 
the  input. 

2.6  PREPALOC  Internal  Common  Blocks 

All  common  blocks  used  internally  by  PREPALOC  are  given  in  table  1.  For 
definition  of  common  blocks  that  communicate  with  the  COP,  see  Program 
Maintenance  Manual,  Volume  I. 
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Table  1.  Module 
(Part 

PREPALOC  Common  Blocks 

1  of  2) 

COMMON 

VARIABLE 

OR  ARRAY 

DESCRIPTION 

CRLENGTH 

CRLENGTH (30) 

Total  length  of  a  penetration  corridor; 
indexed  by  corridor  number 

DISTEF 

DISTEF(50) 

Total  length  of  a  depenetration  corri¬ 
dor;  indexed  by  corridor  number 

DISTEG(50) 

Total  length  of  a  depenetration  corridor 
plus  distance  to  nearest  recovery  base; 
indexed  by  corridor  number 

DREC0V(50,4) 

DESIG  of  recovery  base  associated  with 
depenetration  corridor;  indexed  by 
corridor  number,  up  to  four  bases 

DISTNA(50,4) 

Distance  from  depenetration  corridor  to 
recovery  base;  indexed  by  corridor  num¬ 
ber,  up  to  four  bases 

INDXA(50, 4) 

Ordinal  of  bases  ordered  on  distance 
from  depenetration  corridor;  indexed  by 
corridor  number,  up  to  four  bases 

FIXOFF 

OFFSW 

Switch  to  indicate  if  fix  is  to  be  off¬ 
set  (.true.)  or  not  (.false.) 

GEOREF 

IPNRF 

IDS  reference  code  of  penetration  corri¬ 
dor  header 

IDPRF 

IDS  reference  code  of  depenetration 
corridor  header 

IFXREQ 

IFXREQ(250) 

Fixed  assignment  requests;  indexed  by 
group 

IFXHON(250) 

Fixed  assignment  requests  honored; 
Indexed  by  group 

IFXSW 

F1XSW 

Fixed  assignment  switch:  true  if  there 
are  assignments  in  this  run,  false 
otherwise 

NTBFLG 

NTBFLG (9) 

Alphanumeric  value  to  indicate  source 

of  values  in  general  gaming  parameters. 
May  be  DEFAUL  for  default,  UNCHNG  for 
changed  during  this  run 
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Table  1.  (Part  2  of  2) 


ASSOCIATED 

COMMON 

VARIABLE 

OR  ARRAY 

DESCRIPTION 

RFPOINTS 

RFLAT(20) 

Latitutude  of  a  refuel  point 

RFL0NG(20) 

Longitude  of  a  refuel  point 

MREGN 

Number  of  refuel  points 

SETCOM 

NTBCH 

General  gaming  parameter  switch:  true 
if  changes  were  made,  false  otherwise 

DBCH 

Target  parameter  switch:  true  if  changes 
were  made,  false  otherwise 

NDBCH 

Number  of  target  parameter  change  re¬ 
quests  written  onto  file  25  by  SETRD 

RECALC 

RECALC  mode  switch:  true  if  RECALC  mode 
active,  false  otherwise 

TSTUFF 

Communicates  with  TOFM 

XTOFMIN 

Minimum  time  of  flight 

XCMISS 

Missile  flight  parameter 

XRNGMIN 

Minimum  range 

XRANGE 

Missile  maximum  range 

WEPCOM 

IWRCD(250) 

Weapon  group  record  IDS  reference  code; 
indexed  on  group 

IWTYP(250) 

Weapon  group  type  index 

IWCLS(250) 

Weapon  group  class  indicator  (1-missile, 
2-bomber,  3-salvoed  missile) 

NMWTYP(120) 

Weapon  type  name;  indexed  on  type 

IWMIRV(250) 

MIRV  indicator  (0-nonMIRV,  1-MIRV) 

NASM(2SO) 

Number  of  ASMs  in  group  payload 

NMGRP 

Number  of  groups 
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2.7  Subroutine  ENTMOD 

PURPOSE:  To  control  overall  flow  of  processing 

ENTRY  POINTS:  ENTMOD  (first  subroutine  called  when  overlay  link 

PREP  is  executed) 

FORMAL  PARAMETERS:  None 

COMMON  BLOCKS:  C15,  C30,  ERRCOM,  FIXOFF,  IFXREQ,  IFXSV,  NTBFLG , 

SETCOM 

SUBROUTINES  CALLED:  CHGBAS ,  FIXWEP,  GEOPREP,  HDFND,  INSGET,  MODFY, 

PRNPRP,  RETRV,  SETRD,  WEPPREP 

CALLED  BY:  COP 

Method : 

The  basic  method  employed  by  the  driver  routine  of  PREPALOC  is  to  check 
for  particular  input  clauses  in  a  set  order  and  carry  out  functions 
based  on  their  presence.  As  a  preliminary  step,  the  IDS  error  code 
"R04"  is  set  as  acceptable  so  that  the  retrieval  of  nonexistent  DESIG 
values  will  not  cause  termination  and  the  number  table  (NUMIBL)  is  re¬ 
trieved.  The  input  is  then  scanned  for  a  RECALC  clause,  if  it  exists 
the  general  gaming  parameters  (INITSTRK,  CORMSL,  etc.)  are  set  to  their 
default  values  and  GEOPREP  and  WEPPREP  are  called. 

The  subroutine  now  looks  for  any  SETTING  clauses.  Subroutine  SETRD  is 
called  for  each  occurrence.  If  any  SETTING  clause  called  for  changes 
to  the  target  parameters  (VALUE,  MINKILL,  MAXKILL  or  IDHOB) ,  CHGBAS  is 
also  called.  Then,  the  input  is  scanned  for  FIX  clauses  and  FIXWEP 
called  at  each  occurrence.  Finally,  PRNPRP  is  called  to  produce  any 
standard  or  optional  print. 

Subroutine  ENTMOD  is  illustrated  in  figure  3. 


Figure  3.  (Part  2  of  4) 
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Figure  A.  (part  6  of  9) 
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Figure  4.  (Part  7  of  9) 


2.9  Subroutine  FIXWEP 


PURPOSE: 


ENTRY  POINTS: 


FORMAL  PARAMETERS: 


COMMON  BLOCKS: 


SUBROUTINES  CALLED: 


To  create  fixed  weapon  assignments 
FIXWEP 

ICPT  -  index  of  FIX  clause 

CIO,  C15,  C30,  FIXOFF,  IFXREQ,  OOPS,  TSTUFF, 
WEPCOM,  ZEES 

DIRECT,  DISTF,  HDFND,  HEAD,  INSGET,  MODFY,  NEXTTT, 
RETRY ,  STORE,  TOFM,  WEPIN 


CALLED  BY: 


ENTMOD  (of  PREPALOC) 


Method : 

The  first  step  in  the  FIXWEP  process  is  to  retrieve  the  target  list 
header  (TARNUM)  and  initialize  break-point  tables.  These  break-point 
tables  are  used  to  speed  retrieval  of  records  on  the  LISTXX  chain.  By 
saving  the  direct  reference  code  of  up  to  500  records  evenly  spaced 
along  the  chain,  any  record  may  be  retrieved  by  beginning  the  search 
at  the  nearest  break-point  record  that  precedes  it. 

The  remainder  of  the  process  is  driven  by  the  FIX  clause  which  caused 
the  call.  The  FIX  clause  is  viewed  as  a  series  of  subclauses  each  of 
which  contains  a  subject  and  one  or  more  objects.  The  subject  is  a 
collection  of  attributes  and  each  object  is  a  collection  of  values  for 
those  attributes.  As  a  subject  is  read,  its  component  attributes  are 
used  to  set  up  a  group  of  switches  to  indicate  the  operations  to  be 
carried  out  for  each  object.  The  switches  are: 

ISRNG  -  1  if  one  DESIG,  2  if  two  DESIGs 

IASET  -  1  if  no  set  arrival  time,  2  if  arrival  time  is  set 

FIXWEP  then  reads  in  each  related  object  and  creates  ASSIGN  records  as 
follows.  Each  object  will  set  the  following  values: 

IXGRP  -  the  group  number 

ISTSAL  -  the  selected  salvo  (0  if  non  selected) 

AXRV  -  the  selected  arrival  time  (0  if  IASET  *  1) 

INMA  -  the  number  of  assignments  (set  by  NUMALOC) 

CDSG  -  the  DESIG  of  the  target 

If  there  is  a  DESIG  range  (ISRNG  •  2)  the  values  of  the  range  are  placed 
in  CDSG  in  order  and  the  process  carried  out  for  each. 

For  each  set  of  the  above  values,  FIXWEP  first  checks  to  see  if  the  cur¬ 
rent  group  (IGRP)  is  the  same  as  IXGRP.  If  not,  the  current  group  is 
modified  and  the  new  group  retrieved.  The  target  whose  DESIG  is  CDSG 
is  not  obtained  and  the  corresponding  TARCDE  record  found  through  use 
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of  the  break-point  tables.  Now  the  process  loops  INMA  times.  In  each 
loop  FIXWEP  first  checks  to  see  if  group  resources  are  exceeded.  If 
not,  the  process  proceeds  depending  upon  whether  or  not  the  group  is 
salvoed.  For  a  non-salvoed  group  FIXWEP  simply  stores  a  new  ASSIGN 
record. 

For  a  salvoed  group  the  salvo  is  determined  in  one  of  three  ways.  If 
neither  a  salvo  nor  an  arrival  time  has  been  specified,  a  convenient 
salvo  is  selected.  If  a  salvo  was  specified,  it  is  used.  If  an  arrival 
time  was  specified,  the  time  of  flight  is  calculated  and  the  salvo  set 
to  fulfill  the  desired  time  of  arrival.  In  any  case,  FIXWEP  now  checks 
Co  see  if  a  weapon  is  available  in  the  chosen  salvo.  If  so,  the  ASSIGN 
record  is  stored.  If  not,  FIXWEP  finds  a  salvo  that  is  available  and 
uses  that  salvo  instead  unless  Che  desired  salvo  was  specified.  In  the 
later  case,  FIXWEP  will  search  all  current  assignments  for  Che  weapon 
group  with  the  desired  salvo  and  shift  to  the  available  salvo  that  was 
either  assigned  with  no  specification  or,  failing  to  find  such,  one 
whose  arrival  time  was  specified.  The  new  assignment  may  then  be  made 
to  the  desired  salvo. 

Since  IREG  is  an  attribute  on  both  the  WEPNGP  and  the  TARGET  records 
it  must  be  saved  in  IRGHOLD  after  CALL  DIRECT  ('WEPNGP')  then  restored 
before  CALL  MODFY  ('WEPNGP').  Otherwise,  IREG  may  contain  a  value  from 
the  TARGET  and  degrade  the  IDS  structure  for  the  WEPNGP  record,  since 
IREG  is  a  control  key  for  it. 

Subroutine  FIXWEP  is  illustrated  in  figure  5. 
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Figure  5.  (Part  17  of  23) 


Figure  5.  (Part  18  of  23) 


Figure  6.  (Part  3  of  7) 
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2.11  Subroutine  GEOPREP 


PURPOSE : 

To  calculate  and  store  geographic  data 

ENTRY  POINTS: 

GEOPREP 

FORMAL  PARAMETERS: 

None 

COMMON  BLOCKS: 

CIO,  C15 ,  C30 ,  CRLENGTH ,  DISTEF,  GEOREF,  OOPS 

SUBROUTINES  CALLED: 

DIRECT,  DISTF,  DLETE,  GEOIN,  HDFND,  HEAD,  MODFY , 
NEXTTT,  RETRV,  STORE 

CALLED  BY: 

ENTMOD  (of  PREPALOC) 

Method: 

Depenetration  corridors  are  processed  first.  All  old  TDDIST  and  RDDIST 
records  are  deleted.  Then  a  new  RDDIST  record  is  stored  connecting  the 
corridors  with  their  recovery  bases.  Next,  penetration  corridors  are 
processed.  All  old  TPDIST  records  are  deleted.  Then  all  attributes 
sectors  are  calculated  and  the.  penetration  corridor  record  (PENCRD) 
modified. 

The  subroutine  now  cycles  through  the  target  list.  For  each  target  it 
finds  the  shortest  depenetration  distance,  taking  into  account  the  dis¬ 
tance  to  the  corridor,  the  length  of  the  corridor  and  the  distance  from 
the  corridor  to  the  recovery  base.  A  TDDIST  record  is  created  for  the 
best  depenetration  corridor.  A  TPDIST  record  is  also  stored  connecting 
the  target  list  element  to  every  penetration  corridor. 

Subroutine  GEOPREP  is  illustrated  in  figure  7. 
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Figure  8.  (Part  10  of  12) 
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Figure  10.  Subroutine  WEPIN  (Part  1  of  3) 
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SECTION  3.  ALOC  MODULE 


3.1  Purpose 

The  major  purpose  of  this  module  Is  to  determine  the  optimal  allocation 
of  weapons  to  targets,  using  a  Lagrange  multiplier  technique.  The 
weapons  are  divided  into  weapon  groups  —  each  group  containing  weapons 
of  the  same  characteristics  which  are  geographically  close.  Thus, 
except  for  time  launch  interval  constraints  for  some  "salvoed"  missiles, 
weapons  are  considered  identical  within  groups.  Each  target  is  con¬ 
sidered  individually  for  weapon  assignment.  The  order  of  investigation 
is  the  order  of  the  TARCDE  records  on  the  LISTXX  chain  which  was  deter¬ 
mined  by  the  PLANSET  module.  When  all  targets  have  been  processed, 
another  pass  over  this  chain  begins.  This  process  continues  until  the 
Lagrange  method  has  allocated  all  the  weapons  to  targets.  The  assign¬ 
ments  are  stored  on  a  temporary  file  during  the  process  and  entered  into 
the  Integrated  data  base  or  ASSIGN  records  when  allocation  is  complete. 

The  user  is  able  to  specify  weapon  assignments  through  the  FIX  adverb  to 
the  PREP ALOC  module.  In  this  case  the  ALOC  module  will  optimally  assign 
those  weapons  which  have  not  been  fixed.  In  addition,  there  are  capabil¬ 
ities  which  allow  the  user  to  modify  weapon  range  values,  to  restrict  the 
use  of  MIRV  weapons  by  target  class,  and  to  restrict  the  use  of  any  weapon 
group  by  the  value  of  either  of  the  target  attributes  FLAG  or  CNTRYL. 


3.2  Input 

The  precondition  of  the  integrated  data  base  required  is  that  the  PREFALOC 
module  has  to  have  been  executed.  Furthermore,  there  is  an  optional  input 
file  —  the  Weapon/Target  Data  File.  The  Weapon/Target  Data  file  contains 
the  Information  relating  each  weapon  group  to  each  target.  The  Weapon/ 
Target  Data  file,  if  not  input,  is  created  by  the  FRSTGD  subroutine  on 
pass  one  and  may  be  retained  for  later  executions  of  ALOC.  Two  record 
types  are  alternated  on  the  file.  The  first  record  of  each  pair  contains 
target  data.  The  record  length  is  fixed  and  its  contents  are  shown  in 
table  2.  The  second  record  contains  information  concerning  the  target 
whose  data  are  on  the  previous  record  and  each  of  the  weapon  groups  in 
the  data  base.  The  following  example  will  be  used  to  describe  the  con¬ 
tents  and  format  of  this  record.  Assume  the  target  on  the  first  record 
had  two  hardness  components  (HC1  and  HC2)  and  is  defended  by  terminal 
missiles  defenses.  There  are  five  attacking  groups  as  shown  in  figure 
10.2. 

The  record  contains  four  basic  data  items.  The  first  set  of  data  is  the 
bomber  penetration  corridor  for  each  bomber  group.  This  data  is  packed 
using  five  bits  for  each  entry  as  shown  in  table  2.1.  The  second  set  of 
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Table  2.  Format  of  Target  Data  Record  -  File  Code  15 


Word  Content 


1 

Target  Class 

2 

Target  Type 

3-7 

Time  sensitive  target  input  times 

8-12 

Time  sensitive  target  values 

13 

Country  location 

14 

Country  owner 

15 

Flag 

16 

Number  of  hardness  components 

17 

Number  of  time  components 

,18 

Region 

19 

Complex  indicator  (ICOMPL) 

20 

Latitude 

21 

Longitude 

22 

First  hardness  component  radius  for 

ground  burst 

23 

Second  hardness  component  radius  for  ground  burst 

24 

First  hardness  component  radius  for 

air  burst 

25 

Second  hardness  component  radius  for  ground  burst 

26 

MAXKILL 

27 

MINKILL 

28 

Maximum  value  of  weapons  to  be  used 

on  target 

29 

Index  number 

30 

Target  radius 

31 

SIOP  table  number 

32 

Value  of  target 

33 

Normalized  target  value 

34 

Height  of  burst  indicator 

35 

High  altitude  target  defense 

36 

Low  altitude  target  defense 

37 

Missile  defense 

38 

Number  of  terminal  interceptors 

39 

Target  name 

40 

DESIG 

41 

Value  of  first  hardness  component 

42 

Value  of  second  hardness  component 
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Group 

B1 

Ml 

M2 

B2 

B3 

Type 

Bomber 

Missile 

Missile 

Bomber 

Bomber 

Payload 

Bombs 

MIRV 

MRV 

Bomb  &  ASM 

ASM 

Penetration  Probability 

PB1 

- 

- 

PB2 

PB3 

First  Weapon 

Bomb 

RV 

AY* 

Bomb 

ASM 

Second  Weapon 

- 

- 

RV 

ASM 

- 

*Note:  AY  -  Composite  yield  of  MRV  system  is  used. 


Figure  10.2.  Target-Group  Data 
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Table  2.1.  Format  of  Weapon/Target  Data  Record  -  File  Code  15 
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data  contains  missile  switches.  If  the  missile  is  inactive,  the  switch 
is  set  to  1.  The  missile  data  always  starts  on  a  word  boundary.  The 
third  set  of  data  contains  the  time  of  arrival  for  both  bomber  and  missile 
groups  in  that  order.  The  fourth  set  of  data  consists  of  probabilities 
of  kill.  The  sequence  of  the  data  is  the  probability  of  kill  of  the 
first  weapon  against  the  first  target  hardness  component  PK(1,HC1)  and 
the  probability  of  kill  of  the  first  weapon  against  the  second  target 
hardness  component,  PK(1,HC2).  These  pairs  of  data  are  stored  for  each 
group.  This  information  is  followed  by  the  probability  of  kill  of  the 
second  weapon,-  if  it  exists,  against  each  hardness  component,  PK(2,HC1) 
and  PK(2,HC2).  If  the  target  has  no  terminal  defenses  the  second  group 
of  PKs  is  not  stored  for  missiles.  Data  for  the  second  hardness  com¬ 
ponent  is  also  omitted  if  only  one  hardness  component  of  the  target 
exists.  In  each  of  these  two  cases,  the  data  fields  are  filled  with  the 
next  data  element.  The  record  length  depends  on  the  number  of  bomber  groups, 
number  of  missile  groups,  number  of  target  hardness  components  and  the 
existence  of  missile  defenses.  The  lengths  are  calculated  by  the  following 
expressions. 


LENI(l)  -  LBC+LMS+ (LUND+LS 2U- 1 ) / 2+1 
LENI(2)  -  LBC+LMS+(LUND+2*LS2U-1) / 2+1 
LENI(3)  -  LBC+LMS+ (LDEF+LS 2D- 1) / 2+1 
LENI(4)  -  LBC+LMS+(LDEF+2*LS2D-l)/2+l 
where : 


for  one  hardness  component  and 
no  missile  defense 

for  two  hardness  components  and 
no  missile  defense 

for  one  hardness  component  and 
with  missile  defenses 

for  two  hardness  components  and 
with  missile  defenses 


LBC  ■  (Number  of  bomber  groups  -l)/7+l 
LMS  ■  (Number  of  missile  groups  -l)/36+l 

LUND  *  Number  of  groups  +  number  of  bomber  groups  (+  number  of 
missile  groups  if  there  is  a  random  area  defense) 

LDEF  -  2  times  the  number  of  groups 

LS2U  ■  Number  of  groups  +  number  of  mixed  loaded  bomber  groups 
LS2D  *  LS2U  +  number  of  MRV  missile  groups 


The  Weapon/Target  Data  file  is  created  in  pass  one  and  may  also  be  used 
in  subsequent  runs  (see  RECALC  Mode:  Users  Manual,  UM  9-77,  Volume  III). 
The  format  for  the  Weapon/Target  Data  file  —  file  code  15  —  appears  in 
tables  2  and  2.1. 
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3.3  Output 


As  a  result  of  its  execution  the  ALOC  module  creates  ASSIGN  (record  type 
70)  records  in  the  integrated  data  base.  Further,  the  attribute  NUMALOC 
is  updated  in  every  group  to  reflect  the  actual  number  of  weapons  allo¬ 
cated  from  that  group. 


3.4  Concept  of  Operation 

In  order  to  conserve  storage,  AI0C  is  broken  up  into  two  main  overlays. 
The  first  overlay  is  called  ALCINT.  This  overlay  read  5  any  user  input. 


including  restrictions,  modifications,  print  request  and  so  on.  Further¬ 
more,  the  weapon  data  is  extracted  from  the  integrated  data  base.  The 
second  overlay  (ALCMUL)  controls  the  determination  of  the  allocation. 

The  driver  routine  of  this  overlay  is  MULCON.  Within  the  second  overlay 
there  are  four  segments.  The  first,  FGD,  obtains  target  data  for  pass 
one.  The  second,  SGD,  obtains  target  data  for  passes  two  and  beyond. 

The  third  segment,  STAL,  principally  routines  STALL,  WAD  and  WADOUT, 
allocates  to  targets  without  terminal  ballistic  missile  defenses.  The 
fourth  segment,  DEFAL,  principally  routines  DEFALOC  and  RESVAL  handles 
ballistic  missile  defended  targets. 

3.4.1  Overlay  ALCINT .  The  routines  in  this  overlay  are  straightfcrward 
and  need  little  explanation  beyond  that  below.  However,  the  routine 
DATGRP  places  information  on  a  sequential  file  (file  code  23)  which  is 
used  by  the  FGD  segment.  Table  3  shows  the  format  and  content  of  this 
file. 


3.4.2  Overlay  ALCMUL.  The  design  of  the  weapon-to-target  allocator 
utilizes  a  hierarchy  of  subroutines  operating  at  different  levels  of 
detail.  Figure  11  illustrates  this  hierarchy.  The  major  functions 
associated  with  these  subroutines  are  summarized  below  and  related  to 
the  overall  concept  in  subsequent  paragraphs. 

Subroutine  MULCON  is  the  first  subroutine  in  the  hierarchy  and  is  re¬ 
sponsible  for  the  control  and  adjustment  of  the  Lagrange  multipliers. 
MULCON  monitors  the  rate  at  which  various  classes  and  types  of  weapons 
are  being  allocated  to  the  target  system  and  makes  appropriate  adjust¬ 
ments  in  the  values  of  the  Lagrange  multipliers.  In  this  role,  MULCON 
does  not  need  any  detailed  information  concerning  actual  allocation.  It 
is  concerned  only  with  the  actual  rate  of  allocation  of  the  available 
inventory  as  the  targets  are  processed.  To  obtain  the  assignment  of 
weapons  to  each  successive  target,  MULCON  simply  calls  subroutine  STALL 
(Single  Target  Allocator)  for  targets  without  missile  defenses,  or  sub¬ 
routines  STALL  and  DEFALOC  if  the  target  is  defended.  STALL  and  DEFALOC 
utilize  the  current  values  of  the  multipliers  to  make  an  allocation  to 
the  next  target,  then  return  control  to  MULCON. 

The  data  acquisition  for  the  allocation  process  is  performed  by  the 
FRSTGD  routine  on  pass  one  and  SCNDGD  on  all  other  passes.  Each  of 
these  routines  brings  in  the  proper  records  for  the  next  target 
and  prepares  the  weapon  data  for  that  particular  target.  The  Weapon/ 
Target  Data  File  is  read  (or  on  the  first  pass  in  the  RECALC  mode  cal¬ 
culated).  On  the  first  pass,  FRSTGD  then  creates  a  record  on  file  code 
21  which,  principally,  contains  the  INACTIVE  array  (see  table  4).  If 
the  user  has  requested  range  modification,  FRSTGD  may  also  write  a  re¬ 
cord  onto  file  code  22  in  the  same  format  as  file  code  15  (table  2). 

This  new  file  serves  as  a  source  for  replacement  record  for  file  code 
15.  SCNDGD  reads  each  of  these  files  in  order  to  obtain  the  appropriate 
information. 
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Table  3.  Weapon  Data  File  from  DATGRP 


WORD 

1 

2-4 

5-9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 


DESCRIPTION 

Group  type  index  number 

Logical  flag  restrictions  (99  flags  packed  1  bit  per  flag) 

Logical  country  location  restriction  (150  switches  correspond¬ 
ing  to  codes  in  block  CNCLS ,  packed  1  bit  per  switch) 

Switch  —  true  if  group  is  restricted  MIRV 

Logical  MIRV  restriction  switches  (packed  1  bit  per  switch) 

Range  multiplier 

Refueled  range  multiplier 

Minimum  range  replacement  value 

NALTDLY  for  group 

ALTDLY  for  group 

GLAT  for  group 

GLONG  for  group 

GREFCODE  for  bomber  group,  0  for  missile  group 
GYIELD  for  group 
RANGE  for  group 

CEP  for  bomber  group,  CEP  at  0  range  for  missile  group  (CEP  * 
RANGE  *  (CEP-CEPMIN) / (RANGE-RANGMIN) ) 

SPEED  for  group 

RANGED  for  group 

RANGER  foT  group  (in  DATAMAKE  mode,  PENPROB) 

RNGMIN  for  group 
GREFTIME  for  group 
TOFMIN  for  group 
CMISS  for  group 

Slope  of  CEP  equation  for  missiles,  0  for  bombers 
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Update  Total 
Weight  For 
Next  Target 


Time  To 
Recalculate 
Multipliers  For 
Next  Target? 


No 


*c 


Yes 


Figure  12.  (Part  2  of  2) 
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Figure  12.  Part  XI:  Computation  Loop 
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Again,  after  the  new  allocation  has  been  made,  it  is  added  into  the  run¬ 
ning  sums.  Before  passing  on  to  the  next  target,  the  current  value  of 
the  target  weight  is  revised. 

After  every  two  to  four  targets,  the  Lagrange  multipliers  are  updated 
by  transferring  control  to  the  multiplier  computation  loop  shown  on 
figure  12.  The  error  in  the  rate  of  allocation  for  each  collection  of 
weapons  J  is  estimated.  Three  separate  estimates  are  made  corresponding 
to  differing  rates  of  increase  of  the  target  weights.  If  all  estimates 
have  the  same  sign,  then  a  small  adjustment  of  the  multiplier  in  the 
indicated  direction  is  made. 

The  revised  local  multipliers  are  then  used  to  recalculate  the  Lagrange 
multipliers.  During  the  closing  phase  (PROGRESS  =*  1),  the  local  multi¬ 
pliers  are  not  changed  so  the  actual  multipliers  remain  unchanged.  The 
rate  of  change  of  the  target  weight  is  adjusted  depending  on  the  apparent 
size  of  the  current  error  in  the  allocation  rates. 

Finally  the  progress  of  the  allocation  is  evaluated  and  flags  are  set 
if  the  mode  of  operation  is  to  change. 

3.4. 2. 2  Subroutine  STALL.  STALL  is  basically  a  very  simple  routine. 

It  is  not  responsible  for  computing  payoffs.  The  payoffs  are  computed 
by  WAD  and  sunmarized  for  STALL  by  WADOUT.  Thus,  the  input  to  STALL  con¬ 
sists  only  of  the  summary  data  provided  by  WADOUT.  These  data  consist 
of  only  the  following  variables: 

a.  PPMX  the  maximum  potential  profit  available  if  a  weapon  is 
added,  and  IPPMX  the  index  to  that  weapon 

b.  PVRMX  the  maximum  efficiency  for  any  potential  weapon  that 
could  be  added,  and  IPVRMX  the  index  to  that  weapon 

c.  DFMN  the  minimum  differential  profit  produced  by  any  weapon 
on  the  target,  and  IDPMN  the  index  to  that  weapon. 

Actually,  of  course,  the  profits  and  efficiencies  mentioned  above  are 
based  on  a  modified  payoff  (or  BENEFIT)  computed  by  WADOUT,  which  in¬ 
cludes  the  actual  payoff  from  WAD  together  with  the  premiums  for  stay¬ 
ing  close  to  the  desired  allocation  rates.  Thus,  throughout  the  alloca¬ 
tion,  the  operation  of  STALL  remains  the  same;  it  simply  tries  to  max¬ 
imize  this  modified  payoff.  Changes  in  the  mode  of  the  allocation  are 
thus  accomplished  in  WADOUT  simply  by  changing  the  way  the  payoffs  are 
modified,  so  no  change  in  the  logic  of  STALL  is  required. 
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Figure  13.  Part  I 

(Part  2  of  2) 


Figure  13.  Fart  II:  Fixed  Weapon 
Asaign&ent  Processing 


If  the  IMATCH  parameter  is  used,  VTMIN  and  VTMAX  retain  their  default 
values  until  the  Oth  order  calculation  VTZO  indicates  that  MINKILL  or 
MAXKILL  have  been  reached.  Then  VTMIN  or  VTMAX  is  set  correspondingly 
and  thereafter  operates  as  usual. 

In  any  case,  WADOUT  modifies  its  calculations  so  that  every  weapon  placed 
on  target  destroys  at  least  that  percentage  of  original  target  value  spe¬ 
cified  by  the  user-input  parameter  MINDAMAG. 

3. 4. 2. 5  Subroutine  PREMIUMS.  The  purpose  of  this  subroutine  is  to  com¬ 
pute  the  payoff  premium  required  by  STALL  to  avoid  unduly  large  devia¬ 
tions  from  the  desired  allocation  rate.  During  the  closing  phase,  the 
premium  is  also  used  to  cause  STALL  to  close  in  to  an  allocation  that 
exactly  meets  the  stockpile  constraints.  PREMIUMS  is  called  with  one 
parameter  which  specifies  for  which  group  a  new  evaluation  of  the  pre¬ 
miums  (PREMIUM(G)  and  DPREMIUM(G))  is  desired.  The  call  results  in  the 
replacement  of  the  old  value  of  these  premiums  by  a  new  value. 

3.5  Comnon  Block  Definitions 

The  common  blocks  internal  to  the  ALOC  module  are  shown  in  table  5. 


Table  5.  ALOC  Module  Common  Blocks  (Part  1  of  10) 


BLOCK 

ARRAY  OR  VARIABLE 

DESCRIPTION 

ALERUN 

ALERREST (360,3) 

Estimate  of  allocation  errors 

RUNSUM(360, 3) 

Running  sum  of  target  weight  times 
weapons  allocated 

C33 

NBLN 

Number  of  ballistic  interceptors 

XXYY 

Not  used 

VT 

Target  value  remaining 

TGTWT(3) 

Weighting  values 

PAYOFF 

Target  value  destroyed 

COST 

Sum  of  weapon  values  allocated 

PROFIT 

PAYOFF  -  COST 

DPROFIT 

Change  in  PROFIT  from  last  pass 

WRTEST 

Test  parameter 

CNCLS 

CNTRY(IOO) 

List  of  country  codes 

CLSS (20) 

List  of  target  classes 

NCNT 

Number  of  country  codes 

CONTRO 

WADOP 

WAD  option 

PROGRESS 

Measure  of  allocation  progress 

NPASS 

Allocation  pass 

CORSTF 

COSCOR(30) 

Cosine  of  corridors  orientation  point 
latitude 

DPLAT(30) 

Difference  between  latitude  of  orienta¬ 
tion  and  origin 

DPLONG(30) 

Difference  between  longitude  of  orien¬ 
tation  and  origin 

84 


CH-2 


Table  5.  (Part  1  of  10) 


BLOCK 

ARRAY  OR  VARIABLE 

DESCRIPTION 

CORSTF 

CORLN(30) 

Distance  from  orientation  to  origin 

(cont.) 

• 

TDEFDIST(30) 

Total  precorridor  distance 

TATTRPRE (30) 

Total  precorridor  attrition 

ATTR0C0Z(30) 

Corridor  attrition  rate 

ATTRSUZ(30) 

Corridor  suppression  rate 

HILOATZ(30) 

Ratio  of  high  to  low  attrition 

CRIENGTH(30) 


Corridor  length 


Table  5.  (Part  2  of  10) 


BLOCK 

ARRAY  OR  VARIABLE 

DESCRIPTION 

CORSTF 

0RGLAT(30) 

Origin  latitude 

(cont. ) 

ORGLONG (30) 

Origin  longitude 

DISTCDZ (30) 

Distance  from  origin  to  target 

ATT RAD (30) 

Sum  of  precorridor  and  corridor  attri¬ 
tion 

CR0SSDST(30) 

Perpendicular  distance  from  axis  to 
target 

ENTLAT(30) 

Entry  latitude 

ENTLONG (30) 

Entry  longitude 

DISTDG 

Distance  from  target  to  recovery  base 

MAXCOR 

Number  of  corridors 

CURSUM 

CSALL 

Number  of  targets  assigned  to  all  weap¬ 
ons 

CSREG (5) 

Number  of  targets  assigned  to  weapons 
from  a  given  region 

CSCLAS (2) 

Number  of  targets  assigned  to  weapons 
from  a  given  class 

CSTYPE(IOO) 

Number  of  targets  assigned  to  weapons 
of  a  given  type 

CSGRP(250) 

Number  of  targets  assigned  to  weapons 
from  a  given  group 

CSOTH(2) 

Number  of  targets  assigned  to  weapons 
with  a  given  alert  status 

DEFCOM 

RATM 

Highest  return  rate  from  DEFALOC 

ISALFX(250) 

Storage  for  fixed  salvo  numbers 

NSLC250) 

Number  of  salvoed  weapons  available 

RATE (250) 

Rate  of  return  for  missile  on  defended 
target 

DEFRES 

NOWEP(250) 

Number  of  weapons  assigned  by  DEFALOC 

VTDX 

Surviving  target  value 

NTX(3) 

Terminal  defender  estimates 

FX<3) 

Probability  of  NTX 
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Table  5.  (Part  3  of  10) 


BLOCK 

ARRAY  OR  VARIABLE 

DESCRIPTION 

DYNAMIC 

Contains  allocation 

IG(30) 

Group  assigned 

KORRX(30) 

Corridor  assigned 

RVALX(30) 

Relative  value  of  assignment 

PENX(30) 

Penetration  probability 

TOARR(30) 

Time  of  arrival 

ISAL(30) 

Salvo  number 

NUMFIX 

Number  of  fixed  assignments 

NUM 

Number  of  assignments 

FIRST 

FIRST 

Indicates  first  target  of  pass 

END 

Indicates  end  of  target  list 

F22LSW 

Indicates  file  22  in  use 

FLGSTF 

LFLAG(688) 

Packed  logical  flag  restrictions 

LCNTRY(695) 

Packed  logical  location  restrictions 

NMMRV 

Number  of  restricted  MIRV  groups 

MRNAM(IOO) 

Payload  table  name  of  restricted  MIRV 

LCLAS(64) 

Packed  logical  MIRV  class  restrictions 

RNMUL(250) 

Range  multiplier 

RNRMUL(250) 

Refueled  range  multiplier 

RHMIN(250) 

Range  minimum  replacement 

FORMTT 

INWORD 

Value  which  needs  a  format 

NFORMAT 

Format  for  INWORD 

GRPHDR 

IWGHDR 

IDS  Reference  Code  for  weapon  group 
header 
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Table  5.  (Part  4  of  10) 


BLOCK 

ARRAY  OR  VARIABLE 

DESCRIPTION 

GRPSTF 

Contains  record  from  file  25 

I  TYPE 

Group's  type  index 

LPXAG  (3) 

Flag  restrictions 

LCNTRY (3) 

Location  restrictions 

LMRSW 

Indicate  if  group  is  restricted  MIRV 

LCLASS  (1) 

MIRV  restrictions 

RMUL 

Range  multiplier 

RRMUL 

Refueled  range  multiplier 

RMIN 

Minimum  range  replacement 

GPARAM  (16) 

Contains  the  following  group  param¬ 
eters  in  order:  NALTDLY,  ALTDLY, 
GLAT,  GLONG,  GREFCODE,  GYIELD,  RANGE 
CEP,  SPEED,  RANGED,  RANGER,  RNGMIN , 
GREFTIME,  TOFMIN,  CMISS,  SLOPE 

INITSW 

RECALC 

Indicates  RECALC  mode 

PUNSW 

Indicates  output  of  final  lambdas  is 
desired 

FLAGSW 

Indicates  flag  restrictions 

LOCRSW 

Indicates  location  restrictions 

RMODSW 

Indicates  range  modifications 

MRVRSW 

Indicates  MlftV  restrictions 

PUNIT 

Logical  unit  on  which  final  lambdas 
are  to  be  output 

IRTLEN 

LENI(4) 

Filecode  15  lengths 

KR 

Defense  switch  indicator 

N3TK2 

Switch  which  indicates  if  more  data 
is  needed 
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Table  5.  (Part  4  of  10) 


BLOCK  ARRAY  OR  VARIABLE 

KASG  KGROUP ( 8 ) 

XALAT 

XALONG 

KXORR(15) 

XRVAL(30) 

XPEN(30) 

XARV(30) 

KSALVO(6) 

KHOB(l) 

KFXD(l) 

LACB  LALL 

LAREG(5) 

LACLAS(2) 

LATYPE(IOO) 

LAGRP(250) 

LAOTH(2) 

LGR  LGREF(250) 


DESCRIPTION 

Group 

Target  latitude 
Target  longitude 
Corridor  assignment 
Value  of  target 
Penetration  probability 
Time  of  arrivel 
Salvo  number 

Height  of  burst  indicator 
Fixed  target  indicator 

Lambda  for  all  weapons 
Lambda  for  a  given  region 

Lambda  for  a  given  class 

Lambda  for  a  given  weapon  type 

Lambda  for  a  given  group 

Lambda  for  a  given  alert  status 

Local  Group  Reference  number 


Table  5.  (Part  5  of  10) 


BLOCK 

ARRAY  OR  VARIABLE 

DESCRIPTION 

NALLY 

NALL(250) 

Number  from  group  allocated  on  this 
pass 

RNALL(250) 

Number  from  group  currently  allocated 

NOWPS 

NOALL 

Number  of  weapons  total 

NOREG (5) 

Number  of  weapons  in  a  given  region 

NOCLAS (2) 

Number  of  weapons  in  a  given  class 

NOTYPE(IOO) 

Number  of  weapons  of  a  given  type 

NOGRP(250) 

Number  of  weapons  in  a  given  group 

NOOTH (2) 

Number  of  weapons  with  a  given  alert 
status 

PAYOFF 

OPROFIT 

Profit  from  old  allocation 

SPAYOFF 

Sum  of  all  payoffs 

SUMCOST 

Sum  of  all  costs 

SPROFIT 

Sum  of  all  profits 

PAYSAV 

GSCC(IOO) 

CCREL  for  payload 

GSREL(IOO) 

REL  for  payload 

GSEASM(IOO) 

EXPASM  for  payload 

GS LINT (100) 

LCHINT  for  payload 

NGSWHD(IOO) 

NWHDS  for  payload 

NGSDEC (100) 

DECOYS  for  payload 

GPAYALT (100) 

PAYALT  for  payload 

GYLDASM(IOO) 

ASM  yield  for  payload 

IWHOB(IOO) 

Payload  height  of  burst 

0  ■  for  ground 

1  -  for  air 

2  ■  if  not  preset 
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Table  5.  (Part  6  of  10) 


BLOCK 

ARRAY  OR  VARIABLE 

DESCRIPTION 

PNAV 

GSPKNAV(IOO) 

PKNAV  for  payload 

PREMS 

PREMIUM(250) 

Premium  for  using  weapon 

DPREMIUM(250) 

Premium  for  deleting  weapon 

SUMPREM 

Sum  of  premiums 

TBENEFIT 

Total  benefit 

PREMSV 

PRSV (2) 

Parameter  for  premiums  subroutine 

PRNTCN 

IDO(40) 

Print  request  activation  switch 

INDEXPR(40) 

Print  request  selection  number 

JPASS(40) 

Print  request  first  pass 

JTGTP(40) 

Print  request  first  target 

LPASS(40) 

Print  request  last  pass 

LTGT(40) 

Print  request  last  target 

KTGTFREQ(40) 

Target  print  frequency 

ICOUNT (40) 

Print  request  frequency  counter 

MAXREQ 

Maximum  number  of  requests 

MPRNT 

Number  of  array  entries 

NREQ 

Number  of  requests 

PRTMUL 

DELTEFF 

Increase  in  prof it/VALWPNS 

SDELTEFF 

Sum  of  DELTEFF 

VALWPNS 

Sum  of  all  weapon  values  (lambdas) 

VALERR 

Value  of  surplus  plus  deficit  weapons 

QTKCOM 

QTK 

First  PK 

QTK 

Second  PK 

REFPNT 

RFLAT(IO) 

Refuel  point  latitude 

RFLONG(IO) 

Refuel  point  longitude 

SALVO 

NSALW 

Number  of  salvoed  weapon  groups 

HXSAL(75) 

Maximum  salvo  number  per  weapon  group 
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Table  5.  (Part  6  of  10) 


SALVO 
(cont. ) 


NSALAL(450)  Running  sums  of  salvo  allocation  (six 

words  per  salvoed  group  -  packed  four 
sums  per  word) 

LXIHAVE(50)  Packed  logical  switch  indicating  salvo 

with  weapons 
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Table  5.  (Part  7  of  10) 


BLOCK 

ARRAY  OR  VARIABLE 

DESCRIPTION 

SALVO 

SAVLAM(250)* 

Storage  for  salvoed  weapon  lambdas 

(cont. ) 

(contains  average  payload  difference 
for  bombers  -  AVDE) 

MYSAL(250)* 

Available  salvo  (contains  bomb/ASM 
selection  for  bombers  -  ISETPAY) 

P(250)* 

Balance  parameter  (contains  current 
utilization  of  ASMs  for  bomber  -  FASM) 

ISALPT(250) 

Salvo  index,  indexes  arrays  MXSAL, 
NSALAL  and  LXIHAVE:  0  for  nonsalvoed 
groups 

SMATAD 

SMNOMIRV(3) 

SMAT  parameters  for  non-MIRVs 

SMATMIRV (3) 

SMAT  parameters  for  MIRVs 

SURPW 

SURPWP (250) 

Estimated  weapon  surplus 

TABLE 

TABLE(lOl) 

Table  of  square  root  law  K-factors 

TGTSAV 

TGTLAT 

Target  latitude 

TGTLONG 

Target  longitude 

TGTCLS 

Target  class  name 

VO(2) 

Target  value  per  hardness  component 

WADFIN 

VTP(250) 

Value  remaining  at  target  after  weapon 
added 

DELVT(30) 

Difference  in  surviving  value 

NUMO 

Numbers  of  old  allocations 

IGO(30) 

Group  numbers  of  old  allocation 

*  For  bomber  groups  these  arrays  are  equivalenced  to  arrays  AVDE,  ISETPAY 
and  FASM. 
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Table  5.  (Part  8  of  10) 


BLOCK 

ARRAY  OR  VARIABLE 

DESCRIPTION 

WADFIN 
(cont. ) 

IOP 

Number  of  adds  and  deletes  on  this  tar 
get 

IOPS 

Sum  of  IOP 

WADLOC 

NWP(10) 

Number  of  weapons  in  TOA  set 

VALQ(10) 

Unattritioned  value  at  TOA 

MU(1Q,2) 

Sum  of  means  through  TOA  set 

SIG(10,2) 

Sum  of  variance  through  TOA  set 

V(ll ,2) 

Unattritioned  component  value 

S  (10,2) 

Component  survival  probability  through 
TOA  set 

VS (10,2) 

“(V (N, JH)-V (N+l , JH)  *  S (N, JH) 

VSN (11,2) 

-  VSN(N-l.JH)  +  VS (N-l, JH) 

ITOA (250) 

TOA  index  for  group 

IADDTOA (250) 

1  if  new  TOA  set  required 

SIGP(250,10,2) 

Increase  in  variance  for  TOA  set  if 
weapon  added 

DSIG (250,2) 

Temporary  contribution  of  weapon 

SIGD(30, 10,2) 

Change  in  variance  if  weapon  deleted 

WADOTX 

PVRMX 

Maximum  efficiency 

I PVRMX 

Index  of  weapon  achieving  PVRMX 

PPMX 

Maximum  profit 

IPPMX 

Index  of  weapon  achieving  PPMX 

DPMN 

Minimum  profit 

IDPMN 

Index  of  weapon  achieving  DPMN 

NUMMAX 

Maximum  number  of  weapons  allowed  per 
target 

MV 

Number  on  target 

TPMX 

Largest  potential  profit  . 

NTOA 

Number  of  TOA  sets 

NOTAMAX 

Maximum  TOA  sets 

VTMIN 

Lower  target  destruction  minimum 
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Table  5.  (Part  9  of  10) 


BLOCK 

WADOTX 
(cont . ) 


WADWPN 


WEPSAV 


WPFIX 


ARRAY  OR  VARIABLE  DESCRIPTION 


VTMAX 

ALPHA 

VTEF 

VTZO 

VTO 

STALPRIN 

G 

N 

INACTIVE (250) 

TOA(250) 

TVALTOA (250) 

VTOA(250,2) 

MUP(250,2) 

RISK (6, 250, 2) 

SSIG(250,2) 

MORR(250) 

PEX(250) 

XMUP(250,2) 

I  LAW 

IP(250) 

GSSBL(250) 

IGTYP(250) 

IGLERT(250) 

IGREG(250) 

NWPNS (250) 


Maximum  acceptable  surviving  target 
value 

Factor  on  value  required  to  justify 
VTMAX 

Maximum  of  target  value  remaining  and 
VTMIN 

Total  surviving  target  value 
First  target  component  value 
Stall  print  code 
Group  Number 
Allocation  number 

Active  group  switch  (0  =  Active) 

Time  of  arrival  on  target 

Value  of  target  at  arrival  time 

Value  per  component  at  TOA 

Contribution  of  weapon  to  mean  if  added 

Relative  risk  of  weapon  interaction 

Square  root  of  In  of  SSKP 

Optimal  corridor 

Penetration  probability 

MUP  for  alternate  warhead  -  (ASM  for 
bomber  group,  single  warhead  for  MRV) 

Damage  law  in  use 

Group  payload  index 
Group  SBL 
Group  type  index 
Group  alert  status  index 
Group  region  index 

Number  of  weapons  in  group 
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Table  5.  (Part  10  of  10) 


BLOCK 

ARRAY  OR  VARIABLE 

DESCRIPTION 

WPFIX 

NMTYP (100) 

Type  names 

(cont • ) 

LAH(250) 

Group  lambda 

ITCL(IOO) 

Weapon  type  class  index  (1-missile, 
2-bomber) 

WTS 

WTFAC(3) 

Divide  into  old  weight  for  commeasura- 
bility 

WTRATE(3) 

Rate  of  increase  of  weights 

WTSUM(3) 

Sum  of  weights 

XFPX 

PENALT(30) 

Penetration  probability  for  corridor 
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3.6  Subroutine  ENTMnn 


PURPOSE: 

ENTRY  POINTS: 

FORMAL  PARAMETERS; 
COMMON  BLOCKS: 

SUBROUTINES  CAT.TJrn- 
CALLED  BY: 

Method: 


Entry  module  for  ALOC 

EN1MOD  (first  subroutine  called  when  overlay  ALOC 
is  executed)  3 

None 

C30,  CNCLS,  GRPHDR,  INITSW,  LACB,  NOWPS  PAYSAV 
PRNTCN,  SALVO,  SMATAD,  TABLE,  WEPSAV  WFIX 
OOPS,  IRTLEN  ’ 

INITIAL,  MULCON,  RANSIZ 

MODGET 


±S  called  t0  set  the  record  size  of  file  25  Next  thP 
ALCINT  overlay  is  read  in  and  executed.  If  no  input  error  S  k 2 

pliIri*i.t«ii?S!lJ5:  °Verlay  is  executed.  Finally,  the  final  multi¬ 
pliers  are  written  (or  punched)  if  the  user  has  so  directed. 

Subroutine  ENTMOD  is  illustrated  in  figure  15. 
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Call  PRNPUT 
To  Display 
Input 


. t _ 

Call  DATGRP 
To  Assemble 
Group  Data 


Call  TIMEPRT 
To  Display 
Initialization 
Timer 


RETURN 


Figure  16.  (Part  4  of  4) 
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Figure  18.  (Part  2 

l 


I _ 

Add  Group 
Parameters  To 
File  23  Format 


Figure  18. 


(Part  5  of  6) 
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Figure  18.  (Part  6  of  6) 


3. 8  Subroutine  MULCON 

PURPOSE :  The  primary  purpose  of  MULCON  is  to  adjust  the 

Lagrange  multipliers  during  the  allocation.  It 
also  monitors  the  progress  of  the  allocation,  and 
controls  the  flow  of  the  module  throughout  the 
allocation. 

ENTRY  POINTS:  MULCON 

FORMAL  PARAMETERS:  None 

COMMON  BLOCKS:  ALERUN,  CIO,  C30,  C33,  CONTRO,  CORSTF,  CURSUM, 

DYNAMI,  FIRST,  GRPHDR,  LACB,  NALLY,  NOWS, 

PAYOFF,  PAYSAV,  PREMS,  PRTMUL,  REFPNT,  SALVO, 

SURPW,  TARREF,  TGTSAV ,  WADFIN,  WADOTX,  WADWPN, 
WEPSAV,  WPFIX,  WTS,  LGR 

SUBROUTINES  CALLED;  ADDSAL,  ASGOUT,  BOMPRM,  DEFALOC,  DIRECT, 

FRSTGD ,  LLINK,  MODFY  NEXTTT ,  PRNTALL,  PRNTCON, 
PRNTNOW,  SCNDGD,  STALL,  TIMEME,  ASCSAV 

CALLED  BY:  ENTMOD  (ALOC) 

Method: 

MULCON  proceeds  by  making  multiple  passes  through  the  target  list  until 
the  correct  number  of  weapons  have  been  allocated.  During  this  process, 
the  Lagrange  multipliers  are  modified  to  force  the  allocation  to  converge 
to  the  given  stockpile. 

The  flow  of  operations  in  MULCON  is  illustrated  in  figure  29.  The  first 
sheet  represents  the  flow  in  five  parts  of  the  program;  Part  I,  the  in¬ 
itialization  phase;  Part  II,  first  pass  processing;  Part  III,  the  main 
flow;  Part  IV,  processing  after  allocation;  and  Part  V,  multiplier  ad¬ 
justment.  The  flow  diagrams  are  annotated  with  statement  references  to 
the  FORTRAN  listing  and  are  in  sufficient  detail  to  be  largely  self- 
explanatory.  In  the  following  sections  comments  are  made  only  where  the 
flow  diagrams  require  some  explanation. 

Part  I:  Initialization 

The  routine  begins  by  calling  utility  subroutine  TIMEME  to  initialize 
the  clock.  NPASS  and  ITGT  are  initialized  so  that  the  print  control 
routine  PRNTCON  will  know  what  prints  to  activate,  and  then  PRNTCON  is 
called  to  do  so.  A  call  for  a  summary  print  of  starting  multipliers 
(PRNTALL (11))  follows. 


First  routine  of  overlay  ALCMUL. 
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The  sensitivity  of  the  allocation  errors  to  the  value  of  the  local  multi¬ 
pliers  is  estimated.  With  a  fairly  large  value  for  the  linear  premium 
PRM  all  sensitivities  are  about  the  same  independent  of  the  size  of  the 
group;  thus,  PARTIAL(J)  is  simply  set  to  -1.0. 

RDNSUM  and  WTSUM  are  initialized  as  if  the  starting  pseudoallocation 
were  exact,  but  the  weight  given  to  that  allocation  is  reduced  by  the 
square  root  of  the  number  of  targets  so  that  it  does  not  take  too  long 
for  data  on  actual  allocation  rates  to  produce  a  significant  effect  on 
the  estimated  rates. 

Part  II;  First  Pass  Processing 

This  section  deals  with  the  target  processing  on  the  first  pass.  On  this 
pass  FRSTGD  is  called  which  processes  the  target,  brings  in  the  weapon 
data  and,  if  necessary,  calculates  the  Weapon/Target  Data  File.  FRSTGD 
also  reads  in  any  fixed  allocation.  Finally,  the  pseudoallocation  is 
removed  from  the  running  target  weights. 

Part  III:  Main  Flow  (After  First  Pass) 

This  section  deals  with  the  target  processing  on  passes  after  the  first. 
First  a  series  of  calls  to  PRNTALL  performs  intertarget  prints.  Next 
PRNTCON  is  called  to  reset  print  controls.  Then,  if  this  is  not  the 
first  pass  processing  continues.  SCNDGD  is  called  to  process  target 
data,  bring  in  the  weapon  related  data  and  read  in  the  old  allocation. 

If  PROGRESS  is  2  the  variable  IVERIFY  is  checked.  If  this  is  a  verifi¬ 
cation  pass  the  limits  are  checked  and  if  not  achieved  the  pass  con¬ 
tinues.  If  all  processing  is  complete  PRNTNOW  is  called  for  final  print. 
Then  the  weapon  group  chain  is  cycled  and  the  attribute  NUMALOC  set 
equal  to  RNALL  for  each  group. 

However,  if  processing  is  to  continue,  the  termination  section  is  ignored 
and  the  program  gets  ready  to  generate  a  new  allocation  to  replace  the 
one  just  read.  On  the  first  pass  the  old  allocation  is  a  pseudo-alloca¬ 
tion  and  the  replacement  is  done  elsewhere  (see  Part  II).  The  replace¬ 
ment  is  accomplished  by  removing  the  contribution  of  the  old  allocation* 
to  all  running  sums  before  the  new  allocation  is  generated.  This  can  be 
done  in  this  simple  way  because  the  values  of  COST,  PAYOFF,  PROFIT, 

TGTWT,  and  DPROFIT,  then  in  memory,  are  those  just  read  in  from  the 
TARCDE  record  and  so  they  correspond  to  the  old  allocation.  Since  the 
quantities  RUNSUM  and  WTSUM  were  divided  by  WTFAC  at  the  end  of  the 
last  pass,  the  old  TGTWT  must  also  be  divided  by  this  same  factor  to 
make  it  commensurate  before  it  is  subtracted  out. 

The  reason  that  REVCOST  must  be  computed  is  that  the  values  of  the  multi¬ 
pliers  have  probably  changed  (unless  PROGRESS  ■  1.0)  since  the  prior 


The  fixed  weapons  are  ignored  because  they  do  not  contribute  to  the 
running  sums. 
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Figure  29.  Part  III:  (Part  3  of  3) 
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Flgurt  29.  Part  IV:  (Part  2  of  4) 


170 


Figure  29.  Part  IV:  (Part  4  of  4) 


3.8.1  Subroutine  ADDSAL 

PURPOSE :  This  routine  updates  the  stockpile  for  salvoed 

weapons 

ENTRY  POINTS:  ADDSAL 

FORMAL  PARAMETERS :  IGP  -  group  number 

IOPT  -  option  code 

NUR  -  index  to  1SAL  array 

ISALIN  -  salvo  number 

COMMON  BLOCKS:  DYNAMI ,  SALVO 

SUBROUTINES  CALLED:  None 

CALLED  BY:  MULCON,  STALL,  WAD,  DEFALOC 

Method: 

If  the  weapon  is  a  nonsalvoed  weapon  but  a  bomber,  ISAL  is  set  to  indi¬ 
cate  whether  weapon  is  a  gravity  bomb  or  ASM.  Otherwise  the  variable  IDIFF 
is  set  depending  upon  the  option.  From  this  the  number  in  NSALAL  (packed 
4  per  word)  is  either  incremented  or  decremented. 

Subroutine  ADDSAL  is  illustrated  in  figure  30. 
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Figure  30.  Subroutine  ADDSAL 
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3.8.2  Subroutine  ASGOUT 


PURPOSE : 


To  update  allocation  assignment  records  in  the 
Integrated  data  base 


ENTRY  POINTS: 
FORMAL  PARAMETERS; 
COMMON  BLOCKS: 


ASGOUT,  ASGIN,  ASGSAV 
None 

CIO,  C30,  DYNAMI,  PAYSAV,  TGTSAV,  WEPSAV,  KASG, 
LGR 


SUBROUTINES  CALLED:  DIRECT,  DLETE,  MODFY,  NEXTTT,  STORE,  HEAD,  GLOG, 

SLOG 


CALLED  BY:  MULCON 


Method : 

First  a  logical  switch  is  set  for  each  new  weapon  assignment  to  indicate 
it  is  unassigned.  Next  each  old  assignment  is  compared  to  the  new  assign¬ 
ments  to  see  if  all  values,  except  RVAL,  of  the  old  assignment  are  equal 
to  a  new  assignment.  If  they  are,  the  RVAL  attribute  is  deleted.  Finally 
an  ASSIGN  record  Is  created  for  ail  the  new  assignments  for  which  there 
is  no  match. 

ASGOUT  contains  the  internal  array  LGREF.  This  array  is  initialized 
to  0.  Then,  each  time  an  assignment  is  made  to  a  group,  the  reference 
code  of  the  ASSIGN  record  is  stored  in  LGREF.  Prior,  to  the  storing  of 
any  ASSIGN  record,  if  LGREF  is  non-zero,  the  last  assignment  to  the  sub¬ 
ject  group  which  was  stored  is  retrieved  and  then  the  new  ASSIGN  record 
reference  code  is  stored  in  LGREF.  LGREF  is  set  to  zero  only  if  the 
reference  code  found  in  LGREF  is  the  same  as  that  for  the  record  to  be 
deleted. 

Subroutine  ASGOUT  has  3  entry  points: 

o  ASGOUT:  Writes  assignment  data  to  temporary  file  25  after  each 
pass  during  an  assignment. 

o  ASGIN:  Reads  assignment  data  from  temporary  file  25  before 

each  pass  during  an  assignment. 

o  ASGSAV:  When  assignment  is  complete,  writes  assignment  data  to 
the  IDS  base. 

Subroutine  ASGOUT  is  illustrated  in  figure  31. 
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3.8.3  Subroutine  BOMPRM 


The  purpose  of  this  routine  is  to  maintain  the 
array  containing  the  fraction  of  weapons  all  alio* 
cated  from  each  group  which  are  ASMs. 

BOMPRM 

IDIFF  ■  -1  if  weapons  are  being  deleted 
+1  if  weapons  are  being  added 

C33,  DYNAMI,  HALLY,  PAYSAV,  SALVO,  WEPSAV 

None 

MULCON,  SCNDGD 


This  routine  merely  updates  the  ASM  fraction  array  FASM  in  common  block 
/SALVO/.  The  important  local  variables  are: 

DONE  “  a  logical  array  set  true  if  a  weapon  has  already  been 

processed  to  update  FASM 

TOTW  *  total  number  of  weapons  allocated  from  a  group  on  the 

target 

TASM  “  number  of  ASMs  from  a  group  allocated  on  the  target 

FASM(G)  is  the  fraction  of  currently  allocated  weapons  from  group  G 
which  are  ASMs. 

The  factor  FASM  is  updated  whenever  the  state  of  the  allocation  changes. 
These  changes  occur  when  allocations  from  a  previous  pass  are  removed 
and  when  the  allocation  from  the  present  pass  is  output.  Thus,  BOMPRM 
is  called  from  subroutine  MULCON  on  each  target,  and  by  subroutine  SCNDGD 
for  each  target  after  the  first  pass. 

Subroutine  BOMPRM  has  one  formal  parameter  IDXFF.  If  weapons  are  being 
removed  (previous  pass's  allocation),  then  the  value  of  IDIFF  is  -1.  If 
weapons  are  being  added,  then  IDIFF  is  equal  to  +1.  In  subroutine 
SCNDGD,  the  call  to  BOMPRM  with  IDIFF  equal  to  -1  is  made  after  reading 
the  last  pass  allocation,  just  prior  to  the  update  of  the  running  allo¬ 
cation  suma.  The  call  from  MULCON  with  IDIFF  equal  to  +1  is  made  just 
prior  to  the  running  sum  update. 

Upon  entry  to  subroutine  BGMPSM,  the  routine  checks  variable  NBLN  in 
/C33/.  If  this  variable  is  negative,  the  allocation  in  /DYNAMI/  was 
made  by  subroutine  DEFALOC  and  contains  no  bomber  weapons.  In  this 
case,  the  subroutine  returns  with  no  further  processing.  If  the 


PURPOSE: 

ENTRY  POINT: 

FORMAL  PARAMETERS: 

COMMON  BLOCKS: 
SUBROUTINES  CALLED: 
CALLED  BY: 

Method : 
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allocation  was  made  by  subroutine  STALL,  the  IG  array  of  /DYNAMI/  is 
checked.  For  each  allocation  in  this  array,  the  variable  KORRX  of 
/DYNAMI/  is  checked  to  determine  if  the  weapon  is  a  bomber.  If  not, 
the  next  entry  in  the  IG  array  is  checked.  If  the  weapon  is  a  bomber, 
then  the  value  of  GSEASM  for  the  group  is  checked  (in  /WEPSAV/).  If 
GSEASM  is  equal  to  zero  or  one,  then  FASM  is  set  to  GSEASM  and  process¬ 
ing  continues  with  the  next  entry  in  the  IG  array.  Otherwise,  FASM  is 
updated  for  the  group. 

The  ISAL  array  of  /DYNAMI/  contains  the  indicator  of  bomb  or  ASM  alloca¬ 
tion  (for  bomber  groups  only.  This  array  is  defined  differently  for  mis¬ 
sile  groups).  If  the  value  is  zero,  a  gravity  bomb  was  allocated.  A 
value  of  one  signifies  the  use  of  an  ASM. 

The  total  number  of  weapons  from  the  group  which  are  currently  allocated 
is  kept  in  array  RNAIX  of  common  /NALLY/.  This  array  is  updated  twice 
for  each  target,  just  following  the  call  on  BOMPRM. 

Using  these  variables,  the  value  of  FASM  is  updated  as  follows: 

Define:  TASM  -  number  of  ASMs  allocated  from  group  G  (as  determined 
from  the  ISAL  array)  on  current  target 

TOTW  *  number  of  weapons  allocated  from  group  G  on  current 
target 

Then, 

(FASM  *  RNALL)  +  TASM 

FASM  -  - °  ■— 

new  SMALL  +  TOTW 

Note  that  the  variables  FASM  and  RNALL  in  the  above  equation  are  arrays 
Indexed  by  the  group  number  G. 

Figure  32  displays  the  logic  of  subroutine  BOMPRM. 
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Figure  32.  (Part  2  of  2) 
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3.8.5  Subroutine  PRNTALL 


PURPOSE :  This  routine  provides  a  way  of  calling  the  print 

subroutine  PRNTNOW  that  is  conditional  on  the 
print  control  flags  set  by  PRNTCON. 


ENTRY  POINTS: 

FORMAL  PARAMETERS: 
COMMON  BLOCKS: 
SUBROUTINES  CALLED; 
CALLED  BY: 


PRNTALL 

IOPT  -  Print  option  number 
C30,  CONTRO,  PRNTCN 
PRNTNOW,  TIMEME 

MULCON,  WAD,  WADOUT,  FRSTGD,  RESVAL,  DEFALOC, 
SETPAY 


Method : 

To  provide  convenient  control  over  prints  in  program  ALOC  almost  all 
print  statements  are  contained  in  subroutine  PRNTNOW .  They  are  acti¬ 
vated  by  calling  PRNTNOW (IOPT)  for  the  appropriate  print  option  IOPT. 

If  it  is  desired  to  place  the  print  under  data~input  control  so  that 
the  print  will  not  appear  unless  a  specific  print  request  is  included 
in  the  data  deck,  this  can  be  accomplished  by  calling  PRNTNOW  via  a 
call  on  PRNTALL.  PRNTALL  executes  Che  request  on  PRNTNOW  only  if  the 
print  control  subroutine  PRNTCON  has  set  the  corresponding  print  con¬ 
trol  flag  IDO (IOPT )  active  (i.e.,  ■  3}. 

For  each  call  PRNTALL  first  checks  to  see  if  the  print  has  been  set  ac¬ 
tive  by  PRNTCON.  If  not,  it  immediately  RETURNS  (statement  2)  to  mini¬ 
mize  time  wasted  on  inoperative  print  calls. 

If  the  particular  print  is  active,  PRNTALL  immediately  calls  TIMEME 
(statement  1)  to  stop  the  clock  which  records  active  time  in  the  pro¬ 
gram.  This  makes  it  possible  to  do  a  test  run  with  an  unusual  number 
of  prints  and  still  obtain  a  valid  estimate  of  what  the  running  time 
would  be  without  such  prints.  After  the  call  on  PRNTNOW,  PRNTALL  re¬ 
activates  the  clock  before  returning  to  the  main  program. 

Before  each  print  option  (except  26) ,  PRNTALL  prints  a  heading  identify¬ 
ing  the  optional  print. 

Subroutine  HtNTALL  is  illustrated  in  figure  34. 
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Figure  35.  Subroutine  PRNTCON 
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3.8.7  Subroutine  PRNTNOW 

PURPOSE :  To  produce  optional  printed  output. 

ENTRY  POINTS:  PRNTNOW 

FORMAL  PARAMETERS:  IOPT  -  Print  option  number 

COMMON  BLOCKS:  ALERUN,  C30,  C33,  DYNAMI,  NALLY,  DAYSAV,  PREMS, 

PRTMUL,  SALVO,  SURPW,  TGTSAV,  LACB,  PAYOFF,  WADFIN, 
WADOTX,  WADWPN,  WEPSAV,  WPFIX,  WTS,  NOWPS, 

SUBROUTINES  CALLED:  ABORT,  PRNTOD,  PRNTOF,  PRNTOS,  TIMEME 

CALLED  BY:  PRNTALL,  MULCON 

Method: 

The  formal  parameter  IOPT  determines  which  print  is  produced.  The  re¬ 
sult  of  the  alternatives  appear  in  the  Users  Manual,  UM  9-77,  Volume  III. 
Options  1,  12,  13,  26,  27,  and  28  require  a  subroutine  be  called  which 
contains  the  print  function. 

Subroutine  PRNTNOW  is  illustrated  in  figure  36. 
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SUBROUTINES  CALLED: 


CALLED  BY 


Method: 


This  routine  processes  each  target  in  target  number  order.  Each  call 
causes  the  next  target  to  be  retrieved.  Since  each  record  on  the  tar¬ 
get  list  points  to  either  a  target  or  a  complex  record,  the  next  step 
in  the  process  is  to  retrieve  the  remainder  of  the  target  data.  Next 
the  weapon  data  is  acquired.  This  process  depends  to  a  great  extent 
upon  whether  the  user  has  saved  file  15  from  a  previous  run  of  ALOC. 

If  so,  this  file  is  read  in  and  unpacked.  If  not,  the  file  is  created 
by  cycling  through  the  weapon  groups  and  calculating  the  various  needed 
quantities.  Much  of  the  group  data  needed  for  this  process  is  contained 
on  file  23  where  it  was  stored  by  DATGRP.  If  the  user  has  specified 
range  modifications,  any  information  which  differs  from  that  on  file 
15  is  written  on  file  22  in  the  same  format. 

During  this  process,  the  INACTIVE  array  is  set.  This  array  has  an  entry 
for  each  group  and  is  either  set  to  0  or  100.  0  implies  that  the  group 

is  available  for  allocation  to  the  target.  100  indicates  that  the  group 
is  unavailable  for  one  of  several  reasons:  target  out  of  range,  time 
decay  requirements,  and  flag  location  and  MIRV  restriction.  This  array 
is  written  onto  unit  21. 

The  final  step  is  to  read  in  any  fixed  assignments  to  the  target  and  up¬ 
date  the  assignment  records. 

Subroutine  FRSTGD  is  illustrated  in  figure  38. 
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First  subroutine  of  segment  FGD. 
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Figure  38.  Subroutine  FRSTGD  (Part  1  of  11) 


Figure  38.  (Part  2  of  11) 
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Figure  38.  (Part  3  of  11) 
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Figure  38.  (Part  4  of  11) 
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3.9.1  Subroutine  CRDCAL 

PURPOSE:  Calculate  penetration  probability  and  pick  bomber 

corridor 

ENTRY  POINTS:  CRDCAL 

FORMAL  PARAMETERS:  JORR  -  Corridor  picked 

XEN  -  Penetration  probability 
ZOA  —  Time  of  arrival 
ICLSS  -  Class  index  of  group 

COMMON  BLOCKS:  C30,  CORSTF,  PAYSAV,  REFPNT,  TGTSAV ,  WADWFN, 

WE PS A V ,  XFPX 

SUBROUTINES  CALLED:  DISTF,  EXP,  TOFM 

CALLED  BY:  FRSTGD 

Method: 

First  a  check  is  made  for  match-ups  of  naval  weapons  to  naval  targets. 
Next,  if  the  group  is  a  missile,  its  corridor  is  set  to  zero  and  its 
penetration  probability  and  time  of  arrival  are  calculated. 

For  bomber  groups,  first  the  PKNAV  and  IFENMO  are  checked  to  see  if 
group  is  restricted  to  corridors  2  and  1  respectively.  Otherwise,  each 
corridor  in  turn  is  examined  to  see  if  it  would  provide  the  best  route 
to  the  target.  Corridors  which  would  make  the  distance  too  great  are 
rejected.  For  the  remainder,  the  penetration  probability  is  calculated 
with  the  bomber  using  low  altitude  flight  as  much  as  possible.  The 
highest  penetration  corridor  is  chosen. 

Subroutine  CRDCAL  is  illustrated  in  figure  39. 


Figure  39.  Subroutine  CRDCAL  (Fart  1  of  6) 
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Figure  39.  (Part  3  of  6) 
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3.9.2  Subroutine  FLGCHK 

PURPOSE:  To  check  flag,  location  and  MIRV  restrictions. 


ENTRY  POINTS:  FLGCHK 

FORMAL  PARAMETERS :  FIND  -  True  if  group  not  restricted 

False  if  group  restricted 


J  COMMON  BLOCKS':  C30,  CNCLS,  GRPSTF,  INITSW,  TGTSAV 

SUBROUTINES  CALLED:  None 


CALLED  BY: 


FRSTGD 


Method: 

This  routine  uses  the  logical  arrays  stored  on  file  25  by  DATGRP.  Each 
type  of  restriction  is  checked  for  correlation  between  the  group  and  the 
current  target. 

Subroutine  FLGCHK  is  illustrated  by  figure  40. 
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Figure  40.  (Fart  2  of  3) 


219 


CH-2 


Figure  40.  (Part  3  of  3) 
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3.9.5  Subroutine  PKCALC 


PURPOSE : 

ENTRY  POINTS: 

FORMAL  PARAMETERS: 

COMMON  BLOCKS: 
SUBROUTINES  CALLED: 
CALLED  BY: 


To  calculate  kill  probabilities 
PKCALC 

STA  -  Array  for  normal  PK 

STB  -  Array  for  second  PK 

ICLSS  -  Class  index  of  weapon  group 

C30,  PAYSAV ,  SALVO,  TGTSAV ,  WEPSAV 

None 

FRSTGD 


Method: 

This  routine  calculates  the  kill  probabilities  (PK)  of  the  weapon  group 
against  the  target  using  a  standard  approximation  to  the  circular  cover¬ 
age  function.  Two  PR's  are  calculated  for  each  hardness  component  of 
the  target.  The  first  PK  is  based  upon  the  yield  of  the  group  as  input 
and  the  group's  CEP.  The  second  PK  varies  according  to  the  class  of 
the  weapon.  For  missiles  it  is  the  first  PK  adjusted  by  the  number  of 
warheads  which  make  up  the  group  yield  (for  any  types  but  MRVs  the  first 
and  second  PR's  are  equal).  For  bombers  the  second  PK  reflects  the  PK 
of  any  ASMs  onboard.  Naval  weapon's  have  their  PR's  input  in  the  attri¬ 
bute  GPKNAV.  Average  Destruction  (AVDE)  is  calculated  for  bombers. 

Subroutine  PKCALC  is  illustrated  in  figure  43. 
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Figure  44.  Sub routine  PRNTOF 


3.9.7  Subroutine  RECON 


PURPOSE:  To  reconstruct  data  for  the  WADWPN  block. 

ENTRY  POINTS:  RECON 

FORMAL  PARAMETERS:  ICLSS  -  class  index  of  group,  NOVAL  -  True  if 

group  is  to  be  inactivated.  False  if  group  is  to 
be  activated 

On 

COMMON  BLOCKS:  C30,  C45,  PAYSAV,  SALVO,  SMATAD,  TGTSAV ,  WADWPN, 

WEPSAV 

SUBROUTINES  CALLED:  ALOG,  SETPAY,  TABLEMUP,  VALTAR 

CALLED  BY:  FRSTGD ,  SCNDGD 

Method: 

Subroutine  SETPAY  is  called  to  select  use  of  gravity  bombs  or  ASMs  from 
the  bomber  groups.  If  ASMs  are  selected  probability  STK  2  is  used  in  the 
calculation  of  the  survival  probability  STK  and  the  MUP  array.  The  MUP 
XMUP,  SSIG,  and  RISK  arrays  are  constructed  according  to  the  required 
formulae. 

The  SMAT  array  is  loaded  depending  upon  whether  the  group  is  a  MIRV. 
Subroutine  RECON  is  illustrated  in  figure  45. 
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3. 9. 7.1  Subroutine  RTAPCK 

PURPOSE:  This  routine  handles  the  packing  and  unpacking  of 

file  15  (22  if  RANGMOD  etc.).  In  addition,  it 
handles  the  target  data  I/O  on  file  15. 

ENTRY  POINTS:  RTAPCK,  RTAPIN,  TGTIN,  TGTOUT 

FORMAL  PARAMETERS:  INFL:  WORD  TO  BE  PACKED  OR  UNPACKED 

JDEF  -  ■  target  defended 

*  target  undefended 
10  -  1  ■  unpacks 

2  -  packs 

COMMON  BLOCKS:  C30,  IRTLEN,  WADWPN,  WEPSAV,  PAYSAV,  WPFIX,  TGTSAV, 

QTKCOM 

SUBROUTINES  CALLED:  None 

CALLED  BY:  SCNDGD,  FRSTGD 

Method : 

RTAPCK  has  4  entry  points: 

o  RTAPCK  -  Either  packs  (10*2)  or  unpacks  (10*1)  the  contents  of 
INFL  using  common  blocks  /WADWPN/  and  /QTKCOM/  as  the  source  or 
destination,  respectively,  of  the  data. 

o  RTAPIN  -  Calculates  the  entries  for  common  block  /IRTLEN/ 

o  TGTLIN  -  Reads  in  target  information  from  file  15. 

o  TGTOUT  -  Writes  out  target  information  to  file  15. 

Subroutine  RTAPCK  is  illustrated  in  figure  45.1. 
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Figure  45.1.  RTAPCK  (Entry  RTAPCK) 
(Part  1  of  8) 


240. 


2 


CH-2 


240.3 


CH-2 


♦ 


BRANCH  ON 
ISW 


•  1  (CORRIDOR) 


•  2  (MISSILE  SWITCH) 


-3  (TIME  OF  ARRIVAL) 


•  4  (PER  PROS) 


-  5  (1ST  PK> 


*  B  (2ND  PK) 


Figure  45.1.  (Pert  3  of  8) 


240.4 


CH 


Figure  45.1.  (Part  4  of  8) 


Figure  45.1.  (Pert  6  of  8) 
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Figure  45.1.  (Part  8  of  8) 
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Figure  45.2.  RTAPCK  (Entry  TGTOUT) 
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Figure  45.4.  RTAPCK  (Entry  RTAPIN) 
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3.9.8  Subroutine  SETPAY 


PURPOSE: 

This  routine 
specify  use 

sets  the  bomber  payload  indicators  to 
of  gravity  bombs  or  ASMs. 

ENTRY  POINT: 

SETPAY 

FORMAL  PARAMETERS: 

None 

COMMON  BLOCKS: 

C30,  CONTRO, 
WEPSAV 

DYNAMI,  PAYSAV,  SALVO,  TGTSAV ,  WADWPN, 

SUBROUTINES  CALLED: 

None 

CALLED  BY: 

RECON 

Method: 

This  subroutine  sets  the  bomber  payload  indicators  (array  ISETPAY  in  com¬ 
mon  block  /SALVO/).  The  setting  calculation  considers  the  allocation 
rate  of  ASMs  (array  FASM  in  /SALVO/),  the  actual  fraction  of  ASMs  in  each 
group  array  GSEASM  in /PAYS AV/,  the  damage  difference  between  the  weapons, 
the  average  damage  difference  (array  AVDE  of  /SALVO/),  and  the  state  of 
allocation  progress  (variable  PROGRESS  in  /CONTRO/). 

Five  local  variables  are  of  some  significance  to  the  calculation.  The 
variable  IPREF  is  used  for  temporary  storage  of  the  payload  indicators. 

A  zero  value  specifies  use  of  a  bomb;  a  value  of  one  specifies  use  of  an 
ASM.  The  variable  DEA  contains  the  expected  damage  (ignoring  the  same 
vehicle  planning  factors)  if  a  bomb  were  used.  Local  variable  CONPAY  is 
used  to  weigh  the  difference  in  allocation  rates  relative  to  the  differ¬ 
ence  in  damage  expectancies.  Local  variable  PAYSENS  is  used  to  calculate 
CONPAY. 

The  variable  CONPAY  is  used  to  reflect  the  importance  of  the  allocation 
rate  difference  relative  to  the  damage  difference.  As  CONPAY  decreases, 
the  allocation  rate  difference  increases  in  importance  relative  to  the 
damage  difference  and  vice  versa.  CONPAY  is  defined  as  follows: 

CONPAY  ■  - PAYSENS - 

(CLOSE  +1.5)  *  PROGRESS 

Thus,  PAYSENS  is  merely  a  multiplicative  factor,  set  internally  by  SETPAY 
to  a  value  of  0.1  which  provides  a  base  value  for  CONPAY.  The  denomina¬ 
tor  of  the  right-hand  term  of  the  equation  represents  the  effects  of  allo¬ 
cation  progress.  As  the  allocation  proceeds,  CONPAY  decreases  to  put  more 
weight  on  the  allocation  rate  differences.  Variable  CLOSE  is  a  user  input 
parameter  which  determines  the  size  of  the  closing  forces.  When  PROGRESS 
equals  one,  CLOSE  increases  to  increase  the  closing  forces.  The  table 
below  displays  the  values  of  the  denominator  for  different  values  of 
PROGRESS.  CLOSE  is  set  to  its  default  value,  1.05,  and  CLOSER  (another 
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input  parameter)  is  set  to  its  default  value,  4. 


Table  CONPAY  Denominator  vs.  PROGRESS 


Denominator 


PROGRESS 


1.02 

1.275 

1.913 

2.55 

6.55 

8.55 


.4 

.5 

.75 

1.00  (initial) 

1.00  (after  one  pass) 
1.00  (highest  value  of 
denominator) 


After  PROGRESS  equals  one,  the  value  of  CLOSE  increases  by  an  amount 
CLOSER  for  each  pass  at  PROGRESS  equal  one.  The  maximum  number  of  passes 
with  this  value  of  PROGRESS  is  1.5. 

When  PROGRESS  is  zero  or  in  a  verification  pass,  the  weapon  with  the 
greater  damage  is  selected  for  allocation.  Also,  the  selection  equa¬ 
tions  are  bypassed  if  a  group  has  only  bombs  or  ASMs. 

Figure  46  displays  the  logic  or  subroutine  SETPAY. 
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START 


RETURN 


Figure  46.  Subroutine  SETPAY  (part  1  of  2) 


Figure  46.  (Part  2  of  2) 
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Subroutine  SCNDGD 


* 


PURPOSE: 

ENTRY  POINTS: 


To  read  in  data  on  second  and  later  ALOC  passes 
over  the  target  list. 

SCNDGD 


FORMAL  PARAMETERS:  None 

COMMON  BLOCKS:  CIO,  C30,  FIRST,  TARREF,  WADWPN,  WEPSAV,  WPFIX, 

PNAV,  IRTLEN 


SUBROUTINES  CALLED: 


CALLED  BY: 


BOMPRM,  DIRECT,  FRSPLT,  HEAD,  NEXTTT,  RECON,  TGTIN, 
RTAPCK,  ASGIN 

MULCON 


Method: 

First,  if  this  is  the  first  target  of  the  pass,  units  15  and  21  and  per¬ 
haps  22  are  rewound.  Then  the  next  target  is  read  and  all  its  data 
accessed.  Next  files  21  and  15  are  read  in.  If  file  22  is  in  use  for 
this  target  it  is  read  with  its  contents  replacing  those  of  file  15.  The 
INACTIVE  array  is  now  loaded  and  the  contents  of  file  15  (22)  unpacked 
into  the  WADWPN  block.  Subroutine  ASGIN  is  called  to  read  the  previous 
assignments  and  store  them  in  the  DYNAMI  block.  Finally,  BOMPRM  is 
called  to  remove  the  assignments  effects  from  ASM  totals  and  RECON  is 
called  for  all  groups. 

Subroutine  SCNDGD  is  illustrated  in  figure  47. 


First  subroutine  in  segment  SGD  --  all  other  routines  appear  in  section 
3.9. 
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INACTIVE 

Array 


Figure  47.  (Part  4  of  4) 
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3. 11  Subroutine  STALL 

★ 

PURPOSE:  This  routine  determines  the  sequence  of  weapon 

additions  and  deletions  required  to  achieve  a 
near  optimum  allocation  of  weapons  to  a  target. 

ENTRY  POINTS:  STALL 

FORMAL  PARAMETERS:  None 

COMMON  BLOCKS:  C30,  C33,  CONTRO,  DYNAMI,  SALVO,  SURPW,  WADFIN, 

WADOTX,  WADWPN 

SUBROUTINES  CALLED:  ADD SAL,  INITSAL,  RESTORE,  WAD 

CALLED  BY:  MULCON 

Method : 

STALL  controls  WAD's  addition  and  deletion  of  weapons  by  setting  the 
values  of  the  following  three  variables: 

WADOP  in  /CONTRO/ 

G  in  /WADOTX/ 

NW  in  /WADOTX/ 

WADOP  has  the  following  options: 

WADOP  ■  1  Initialize  allocation 

WADOP  ”  2  Finalize  allocation  (optional  print  of  results) 

WADOP  *  3  Add  weapon  from  group  (G) 

WADOP  *  4  Delete  (NW)th  weapon  now  on  target 

To  facilitate  monitoring  the  operation  of  STALL,  the  variable  STALPRIN 
is  also  set  to  provide  a  unique  indicator  of  the  position  in  STALL 
where  WAD  is  called.  This  variable  is  printed  under  the  print  option  22. 

The  input  data  for  STALL  consists  of  the  following  six  variables  supplied 
by  WADOUT  in  /WADOTX/ : 

PPMX  and  IPPMX  -  the  maximum  potential  increase  in  effective 
profit  for  any  weapon,  and  the  index  G  to  that  weapon  group, 
respectively. 

*  „ 

First  subroutine  of  segment  STAL. 
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PVRMX  and  IPVRMX  -  the  maximum  effective  efficiency  for  any  weap¬ 
on,  and  the  index  to  that  weapon  group,  respectively. 

DPMN  and  IDPMN  -  the  minimum  marginal  effective  profit  for  any 
weapon  now  assigned,  and  the  index  to  that  weapon  in  the  list  of 
weapons  assigned,  respectively. 

The  flowchart  for  STALL  is  in  four  parts.  The  first  part  contains  the 
setup  and  single  weapon  allocation  phase.  This  phase  provides  a  prompt 
exit  from  STALL  if  the  indicated  allocation  consists  of  one  weapon  or 
less.  This  part  also  includes  a  dummy  version  of  STALL  which  is  used  to 
reproduce  the  prior  allocation  independent  of  current  payoff  data.  This 
option  is  used  in  place  of  the  usual  verification  phase  (when  PROGRESS  * 

2)  if  IVERFY  =  2.  This  mode  of  operation  checks  the  effects  of  an  alter¬ 
nate  level  of  interweapon  correlation,  C0RR2. 

Before  going  into  the  normal  allocation  phase,  the  initial  value  of  the 
time-of-arrival  error  allowance  DELTVAL  is  saved,  so  that  it  can  be  re¬ 
stored  if  it  is  necessary  to  change  it.  This  quantity  determines  the 
maximum  fractional  difference  in  target  value  at  the  time-of-arrival  of 
weapons  that  are  allowed  to  use  the  same  time-of-arrival  bin  in  the  cal¬ 
culations  by  WAD.  If  the  indicated  allocation  would  result  in  an  over¬ 
flow  of  the  available  time-of-arrival  bins,  this  quantity  is  increased 
by  STALL  and  the  allocation  is  reinitialized  for  another  attempt. 

Before  performing  any  operations,  STALL  first  calls  subroutine  INITSAL. 
This  routine  initializes  the  arrays  in  common  block  /SALVO/.  At  the  end 
of  the  routine,  subroutine  RESTORE  restores  the  multipliers  for  the  sal- 
voed  groups. 

The  second  part  of  STALL  processes  the  fixed  assignment  data.  It  puts 
the  weapons  down  on  target.  On  the  first  pass,  the  fixed  assignments 
were  placed  in  array  IG  by  FRSTGD.  In  later  passes,  they  are  in  the 
IGO  array.  The  statements  after  126  determine  the  number  of  weapons  the 
assignment  represents.  If  DEFALOC  has  made  an  allocation  on  the  previous 
pass,  the  number  of  missiles  allocated  from  each  group  is  shown  as  a 
negative  number  in  the  KORRX  array.  If  the  KORRX  entry  is  positive,  there 
is  only  one  weapon  assigned  from  the  group.  STALL  then  checks  the 
INACTIVE  flag  to  see  if  the  weapon  can  reach  the  target.  If  not,  an 
error  message  is  printed  and  the  assignment  request  is  ignored.  If 
statement  443,  WAD  is  called  to  actually  put  the  weapon  on  target.  If 
the  weapon  Is  a  salvoed  missile,  subroutine  ADDSAL  is  called  to  modify 
the  salvoed  weapon  stockpile. 

For  the  fixed  weapons,  no  change  to  the  variable  SURPWP  is  made  as  the 
weapon  is  allocated,  since  this  variable  controls  the  allocation  of  only 
those  weapons  used  by  the  mathematical  allocator. 

If  the  fixed  weapon  cannot  be  allocated  because  of  its  active  flag 
(INACTIVE  (G))  prevents  allocation,  the  error  message  notifying  the  user 
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of  this  fact  lists  the  reason  for  the  inactivity.  The  reason  is  contained 
in  array  MORRX  for  all  inactive  groups. 

The  third  part  of  the  subroutine  provides  an  initial  laydown  of  weapons 
if  multiple  weapons  are  indicated  against  the  target.  As  this  laydown 
progresses  two  types  of  array  overflows  could  occur.  The  overflow  of 
available  time-of-arrival  bins  has  just  been  discussed.  It  results  in 
simply  restarting  the  allocation.  The  other  possibility  is  that  the  total 
number  of  weapons  assigned  could  exceed  the  maximum  number  (30)  permitted. 
In  this  event,  control  is  passed  to  the  refinement  loop  just  as  it  would 
be  in  the  normal  exit  at  the  bottom  of  Part  III. 

Part  IV  independently  checks  for  such  a  potential  overflow  (near  state¬ 
ment  59)  and  if  it  is  threatened,  a  cycle  of  operation  is  generated 
through  the  connector  that  results  in  removal  of  the  least  profit¬ 
able  weapons  (statement  52)  and  replacement  by  the  most  profitable  avail¬ 
able  weapons  (statement  56).  This  sequence  is  terminated  either  by  en¬ 
tering  the  normal  refinement  process  when  the  residual  target  value  is 
reduced  so  that  no  potential  weapons  remain  profitable  (statement  54),  or 
by  using  exit  after  statement  59  if  no  combination  of  30  weapons  can 
be  found  which  will  reduce  the  target  value  sufficiently. 

The  operation  of  the  normal  refinement  loop,  which  cycles  through  the 
branch  at  statement  71  until  the  tests  at  statement  66  are  satisfied,  is 
discussed  in  detail  in  appendix  A  and  will  not  be  repeated  here. 

Figure  48  illustrates  segment  STALL. 
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Pert  II:  Fixed  Weapon 
Assignment  Processing 
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3.11.1  Subroutine  FORMATS 


PURPOSE: 

ENTRY  POINTS: 

FORMAL  PARAMETERS: 
COMMON  BLOCKS: 
SUBROUTINES  CALLED: 
CALLED  BY: 


This  subroutine  determines  the  best  10-column  BCD 
format  for  a  variable. 

FORMATS 

None 

FORMTT 

None 

PRNTOS 


Method: 

The  variable  to  be  formatted  is  INWORD,  the  first  word  in  common  /FORMTT/. 
(This  word  is  equivalenced  to  WORDIN.)  The  best  10-column  format  is  re¬ 
turned  in  NFORMAT,  the  second  word  of  common  /FORMAT/.  The  resulting 
value  of  NFORMAT  can  be  used  in  a  FORTRAN  output  statement  such  as: 

PRINT  NFORMAT,  INWORD. 

The  names  WORDIN  and  INWORD  are  equivalenced  to  allow  correct  specifica¬ 
tion  (real  or  integer)  for  any  possible  input.  The  name  NFORMT  is  in¬ 
ternally  equivalenced  to  the  variable  name  NFORMAT  for  convenience.  In¬ 
ternal  to  FORMATS,  the  absolute  value  of  the  input  variable  is  kept  in 
INABS,  equivalenced  to  ABSIN  for  type  specification  purposes. 

To  determine  if  a  number  is  fixed  or  floating  point,  the  first  12  bits 
of  the  absolute  value  are  tested.  If  a  bit  is  set,  the  number  is  assumed 
to  be  floating  point,  since  all  normalized  floating  point  numbers  have  at 
least  one  bit  set  in  this  range.  Thus,  if  an  integer  quantity  greater 
than  16,777,216  is  input,  it  will  be  treated  as  if  it  were  a  floating 
point  variable.  However,  no  variable  input  from  PRNTNOW  can  have  a  value 
of  that  magnitude,  so  this  restriction  is  never  a  handicap.  Table  6 
presents  the  returned  formats  for  each  range  of  input  values. 

Subroutine  FORMATS  is  illustrated  in  figure  49. 
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Table  6.  Calculated  Formats  for  Variables 


range  of  absolute 

VALUE  OF  VARIABLE,  X 

FORMAT  OF  NEGATIVE 

FORMAT  OF  POSITIVE 

0  or  Integer  Variable 

110 

110 

0  <  X  <  0.0001 

E10.1 

E10.2 

0.0001  *  X  «  0.999999 

F10.5 

F10.5 

0.999999  <  X  s  99.9999 

F10.4 

F10.4 

99.9999  <  X  £  9999.99 

F10.3 

F10.3 

9999.99  <Xs  999999.9 

F10.2 

F10.2 

999999.9  <  X 


F10.1 


F10.2 


3-11.3  Function  LAMGET 


PURPOSE: 


ENTRY  POINTS: 
FORMAT  PARAMETERS: 


COMMON  BLOCKS: 
SUBROUTINES  CALLED: 
CALLED  BY: 

Method: 


»uiPrLaerv^01rL1f„u."5l^s“;^s""«c,''  U8™8e 

LAMGET 

LAM  -  Initial  multiplier 
*  ”  Salvo  balance  variable 

laAL  -  Salvo  number 


DEFALOC,  SA1VAL  (entry  INITSAL) 


~  ££  Its  co„,ct  „3,ee 

gram.  p  citlcation  be  present  in  the  calling  pro- 

Slot’s  °Anrsc  -Itlpli.r 

number  (ISAL).  variable  maintained  by  PUPDT(P),  and  the  salvo 


The  returned  value  LAMGET  is  computed  as  follows: 

LAMGET  -  LAM  -P  *  (ISAL  -1)  *  LAM 
Function  LAMGET  is  illustrated  in  figure  51. 
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Figure  51.  Function  LAMGET 
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3.11.4  Subroutine  PREMIUMS 

PURPOSE:  This  routine  calculates  the  premium  used  by  WADOUT 

in  evaluating  the  benefit  of  using  or  not  using 
weapons  from  specific  groups. 

ENTRY  POINTS:  PREMIUMS 

FORMAL  PARAMETERS :  G  -  An  integer  group  number 

COMMON  BLOCKS:  C30,  CONTRO,  PREMS,  SURPW,  WPFIX 

SUBROUTINES  CALLED:  None 

CALLED  BY:  WAD,  DEFALOC,  MULCON 

Method: 

The  formal  parameter  G  specifies  for  which  weapon  group  the  premium  is 
to  be  recomputed.  Three  modes  are  provided  for  the  computation  of  the 
premiums,  depending  on  the  value  of  PROGRESS. 

PROGRESS  <1.0  A  normal  linear  premium  is  computed  which 
keeps  the  allocator  from  producing  alloca¬ 
tions  with  unnecessarily  large  deviations 
from  the  desired  allocation  rate  (statement 
12). 

PROGRESS  ■  1.0  The  normal  premium  is  augmented  by  a  step 

function  which  strongly  motivates  the  allo¬ 
cator  to  exactly  match  the  stockpile  state¬ 
ment  1). 

PROGRESS  >1.0  The  subroutine  exists  with  a  zero  premium  for 
use  in  verification  allocations  (statement  10). 

Since  the  calculation  of  the  step  function  premium  requires  the  quantity 
SMALLAM  *  .5  *  LAMEF(G)  for  the  lowest  value  of  LAMEF,  this  quantity  is 
evaluated  only  on  the  first  call  of  PREMIUMS  after  PROGRESS  is  1.0  (state¬ 
ment  18).  Thereafter  (since  the  values  of  LAMEF  are  frozen  while 
PROGRESS  -  1.0)  there  is  no  need  to  recompute  this  quantity. 

PREMIUM  is  subtracted  from  the  cost  of  adding  a  weapon.  DPRQilUM  is 
subtracted  from  the  cost  of  deleting  a  weapon. 

Where  a  weapon  surplus  exists  (i.e.,  SURPWP(G) > 0) ,  PREMIUM (G)  >  0.0  and 
DPREMIUM(G) < 0.0.  Where  a  weapon  deficiency  exists  (i.e.,  SURPWP(G) <  0) , 
the  reverse  is  true. 

Figure  52  illustrates  subroutine  PREMIUMS. 
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RETURN 


Figure  53.  Subroutine  PRNTOS 
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3.11.6  Subroutine  SALVAL 


PURPOSE : 


ENTRY  POINTS : 


This  routine  selects  the  preferred  salvo  for  each 
salvoed  missile  group,  saves,  and  restores  the 
appropriate  Lagrange  multipliers. 

SALVAL,  INITSAL,  NEWSAL,  RESTORE 


FORMAL  PARAMETERS:  None 


COMMON  BLOCKS:  C30,  CONTRO,  PAYSAV,  SALVO,  TGTSAV,  WADOTX,  WADWPN, 

WEPSAV,  WPFIX 


SUBROUTINES  CALLED: 
CALLED  BY: 


GLOG,  LAMGET,  VALTAR 
STALL,  DEFALOC,  WAD 


Method : 

There  are  three  entry  points,  INITSAL,  RESTORE,  and  NEWSAL.  Entry 
INITSAL  is  called  at  the  beginning  of  automatic  weapon  allocation  in 
STALL  or  DEFALOC.  First,  it  saves  the  multipliers  LAM  for  each  sal¬ 
voed  group  in  array  SAVLAM.  Then,  for  each  group  and  each  salvo, 

INITSAL  determines  the  salvo  number  and  cost  of  the  best  salvo.  The 
profit  of  each  salvo  is  defined  as 

PRi  -  VTAR1  -TLAMi 

where  VTAR^  is  the  value  of  the  target  at  arrival  time  of  salvo  i  and 
TLAMj_  is  the  value  of  the  multiplier  for  salvo  i  as  determined  by  func¬ 
tion  LAMGET.  The  salvo  with  the  highest  PRj  is  selected  as  the  best 
salvo.  The  salvo  number  is  entered  in  array  MYSAL  in  comnon  block 
/SALVO/  and  the  cost  L^  replaces  the  multiplier  for  the  group  LAM. 

Entry  RESTORE  merely  restores  the  original  values  of  the  multipliers 
(SAVLAM)  back  to  the  multiplier  array  LAM.  Entry  NEWSAL  checks  the 
current  allocation  and  running  sum  for  the  weapon  just  added  or  deleted. 
If  the  salvo  has  been  completely  allocated,  NEWSAL  flags  the  salvo  as 
unavailable.  Entry  SALVAL  is  never  used. 

Entry  INITSAL 

This  entry  is  the  most  complex  of  SALVAL.  The  local  variable  IAM  is  set 
to  "INITSL"  to  flag  the  exit  points  after  statement  420  and  at  statement 
500.  The  first  processing  (DO  loop  to  statement  10)  saves  the  Lagrange 
multipliers  (LAM)  in  the  SAVLAM  array  in  /SALVO/.  The  major  processing 
in  INITSAL  occurs  in  the  DO  loop  to  statement  100  over  all  the  salvoed 
groups.  Within  this  loop,  the  DO  loop  to  statement  300  investigates 
each  salvo  to  determine  its  worth.  Array  IHAVE  in  /SALVO/  is  used  to 
exclude  salvos  from  consideration.  If  IHAVE  (I,  J)  is  false,  then 
aalvo  I  in  group  J  does  not  have  any  available  weapons  and  is  ignored. 

If  weapons  were  available,  a  jump  is  made  to  statement  420  in  entry 
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NEWSAL.  That  section  of  code,  described  later,  determines  if  the  salvo 
to  be  considered  is  overallocated  beyond  its  limit.  If  so,  the  salvo 
is  not  further  considered.  If  the  salvo  is  available,  function  VALTAR 
is  used  to  obtain  VTAR  the  target  value  for  the  salvo.  Function  LAMGET 
is  used  to  calculate  TLAM,  the  multiplier  for  the  salvo.  The  best  salvo 
is  the  one  with  the  highest  positive  difference  between  VTAR  and  TLAM. 
When  this  maximum  is  selected,  MVSAL  in  /SALVO/  is  set  to  the  best  salvo 
number.  LAM  in  /WPFIX/  is  set  to  the  appropriate  multiplier.  The  arrays 
TOA,  TVALTOA,  and  VTOA  (in  /WADWPN/)  are  set  to  correspond  to  the  arrival 
time  of  the  best  salvo. 

Entry  INITSAL  is  illustrated  in  part  1  of  figure  54. 

Entry  RESTORE 

This  simple  entry  point  uses  a  DO  loop  to  statement  1000  to  restore  the 
original  values  of  the  Lagrange  multipliers  (SAVLAM)  to  the  multiplier 
array  (LAM) . 

This  entry  is  illustrated  in  part  3  of  figure  54. 

Entry  NEWSAL 

This  entry  is  used  after  each  allocation  of  salvoed  weapons  to  determine 
if  the  salvo  is  still  available.  Each  salvo  has  a  maximum  limit  of 
overallocation.  Before  PROGRESS  equals  one,  this  limit  is  +225  which 
is  the  largest  number  which  can  be  stored  in  the  NSALAL  array.  (A  de¬ 
scription  of  the  structure  of  the  NSALAL  array  is  contained  in  the  Method 
section  of  subroutine  ADDSAL,  in  this  chapter.)  When  PROGRESS  equals 
one,  a  more  severe  limit  is  imposed  in  order  to  accelerate  closure  to 
the  stockpile.  In  this  case  the  limit  is  zero.  That  is,  a  salvo  is 
available  only  if  it  is  undera 1 located . 

To  exit  from  this  entry,  local  variable  IAM  is  checked.  If  it  is 
"INITSL,"  then  the  original  call  was  to  INITSAL  and  control  passes  to 
statement  220  in  entry  INITSAL.  If  it  is  "NEWSAL"  the  routine  exits. 

Note  that  if  a  salvo  is  unavailable  because  of  limit  on  a  call  to  NEWSAL, 
the  routine  attempts  to  find  the  closest  lower  available  salvo.  If  none 
can  be  found,  then  MTSAL  is  set  negative  to  flag  that  no  salvo  is  avail¬ 
able  for  this  group. 

Entry  NEWSAL  is  illustrated  in  part  4  of  figure  54. 


In  line  function  IHAVE  extracts  information  from  logical  array  LXIHAVE. 
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3.11.8  Subroutine  WAD 

PURPOSE :  This  routine  carries  out  the  addition  and  deletion 

of  weapons  as  specified  by  STALL.  After  each 
change  in  the  allocation  to  a  target  WAD  computes 
the  surviving  target  value  VT ;  for  each  potential 
weapon  G,  the  potential  surviving  value  VTP(G)  if 
a  weapon  for  that  group  were  added;  and  for  each 
weapon  currently  on  the  target,  the  potential  sur¬ 
viving  value  VTD(NW)  if  a  weapon  from  that  group 
were  deleted. 

ENTRY  POINTS:  WAD 

FORMAL  PARAMETERS:  None 

COMMON  BLOCKS:  C30,  C33,  CONTRO,  DYNAMI,  PREMS,  SALVO,  SURPW, 

TGTSAV,  WADFIN,  WADLOC,  WADOTX,  WADWPN,  WEPSAV, 

WPFIX 

SUBROUTINES  CALLED:  ADDSAL,  FMUP,  NEWSAL,  PREMIUMS,  PENT ALL, 

WADOUT 

CALLED  BY:  STALL 

Method : 

The  surviving  target  value  VT  is  given  by: 

J=M  N=NN 

VT  -  Y,  L  [  V  (N,  J)  -V  (N  +  1,J)]  *  S(N,J) 

J=1  N=0 

where 

S(N,J)  =  FMUP  [(MU(N,J)  **  2) / (MU(N, J)  +  SIG(N,J))] 

The  function  FMUP(X)  is  defined  as  follows: 

Exponential  damage  law; 

S (N, J)  *  FMUP(X)  -  exp  (-X) . 

Square  root  damage  law: 

S(N, J)  -  FMUP(X)  -  (1  +>/T")  (exp  (-'/x-)). 

The  index  J  is  over  the  hardness  components,  and  the  index  N  is  over  the 
time  of  arrival  bins.  V(N,J)  is  the  unattrited  value  of  the  Jth  hard¬ 
ness  component  at  the  time  corresponding  to  the  Nth  time-of-arrival  bin. 


However,  the  value  of  VS  would  be  changed  since  the  value  of  V (N  +  1,J) 
should  now  reflect  the  target  value  for  the  new  column  3  and  would  be 
higher  than  the  32.0  now  shown  for  column  3,JH  *  1.  The  value  of  MU 
and  SIG  in  the  new  column  3  would  be  the  same  as  that  in  column  2  ex¬ 
cept  that  MU  would  be  augmented  by  MUP(G.J)  for  the  new  weapon  and  SIG 
would  be  augmented  by  the  cross  term  between  the  new  weapon  and  all  other 
weapons  on  the  target  at  that  time.  The  same  rule  of  course  applies  for 
all  succeeding  time-of-arrival  bins.  In  each  following  column  (includ¬ 
ing  old  column  3,  now  4)  the  previous  value  of  MU  is  increased  by  MUP 
for  the  weapon  added,  and  the  value  of  SIG  is  increased  by  the  cross 
terms  between  all  weapons  previously  on  the  target  and  the  weapon  added. 

Of  course,  the  new  weapon  might  fit  into  one  of  the  existing  time-of- 
arrival  bins.  In  this  case  it  would  be  unnecessary  to  make  a  new  column, 
and  it  would  be  unnecessary  to  recompute  VS  for  the  previous  column.  The 
value  of  MU  and  SIG  would  simply  be  augmented  in  the  corresponding  col¬ 
umn  and  in  all  succeeding  columns  as  before.  Naturally,  after  the  values 
of  MU  and  SIG  are  revised,  the  value  of  S  must  be  recomputed  and  the  VS 
and  VSN  must  be  revised  in  the  columns  affected. 

We  now  recall  that  WAD  is  required  to  provide  potential  weapon  added 
target  value  VTP(G)  for  every  weapon  group  each  time  a  weapon  is  added 
or  deleted.  Obviously  the  calculation  every  time  of  all  the  above  cross 
terms  could  be  very  time  consuming.  Moreover,  precalculation  of  all  in¬ 
dividual  cross  terms  for  250  weapon  groups  would  be  equally  impractical. 

The  technique  adopted,  therefore,  was  to  calculate  cross  terms  for  each 
weapon  group,  but  only  with  the  weapons  already  on  the  target.  Since 
the  process  of  augmenting  the  values  of  SIG  must  be  carried  out  indi¬ 
vidually  for  each  time-of-arrival  bin,  the  resulting  data  are  stored 
by  time-of-arrival  bins.  That  is,  for  each  weapon  group  G  and  each 
time-of-arrival  column,  N,  data  are  stored  which  indicates  the  amount 
by  which  SIG  would  be  increased  in  that  column  if  the  weapon  G  were 
added.  These  data  are  stored  in  the  array  SIGP(G,J,N).  Using  these 
data,  the  effective  augmentation  of  SIG  to  calculate  VTP(G)  for  each 
weapon  group  G  can  be  accomplished  simply  by  adding  SIGP  for  the  appro¬ 
priate  column,  and  the  augmentation  of  MU  is  accomplished  simply  be  add¬ 
ing  MUP  for  each  weapon  G. 

Table  8-A  may  help  to  visualize  how  these  data  are  used.  Each  potential 
group  G  is  tagged  with  the  index  ITQA(G)  of  the  time  of  arrival  column 
it  would  occupy  if  it  were  added.  In  addition,  it  is  also  noted  whether 
the  weapon  would  generate  a  new  column  (IADDTQA~1)  or  share  the  column 
with  the  weapons  already  there  (IADDTQA-0) .  The  situation  illustrated 
here  in  table  8  corresponds  to  the  same  one  illustrated  in  table  7.  No¬ 
tice  that  for  each  weapon  group  G  the  array  SIGP(G,J,NI)  contains  a  sig- 
^  nlflcant  (usually  nonzero)  data  for  NI>IT0A(G)  -LADDTOA(G) .  The  extra  term 
in  the  column  NI  -  IT0A(G)  -1  for  rows  1  and  5,  where  the  weapon  group 
would  require  a  new  column  (IADDTQA*1) ,  is  to  provide  a  term  required 
for  the  new  columns  if  this  weapon  were  added.  Since  addition  of  this 
weapon  would  move  all  columns  (including  its  own)  one  position  to  the 
right,  the  resulting  term  would  be  in  the  proper  position  after  moving, 
even  though  it  is  in  an  incorrect  column  at  present. 
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Table  8.  Illustrating  Quantities  Calculated  for 
Potential  Weapon  Added  and  Deleted 
Payoffs 

A.  Data  on  All  Potential  Weapons 


G 

VTP 

ITOA 

IADDTQA 

J 

DSIG.G.3 

3 

M 

II 

♦— » 

NI=2 

NI=3  NI=4 

(G) 

(G) 

(G) 

1 

2 

0 

1 

.021 

0.0 

.042 

.063 

2 

.103 

0.0 

.206 

.659 

2 

3 

0 

1 

.052 

0.0 

0.0 

.45 

2 

.660 

0.0 

0.0 

4.23 

3 

2 

0 

1 

.274 

0.0 

.548 

1.653 

2 

.780 

0.0 

1.560 

1.880 

4 

3 

0 

1 

.005 

0.0 

0.0 

.010 

2 

.020 

0.0 

0.0 

.040 

5 

4 

1 

1 

.003 

0.0 

0.0 

.007 

2 

.009 

0.0 

0.0 

.018 

B.  Data  on  Weapons  Now  on  Target 

(As  Candidates  for  Possible  Deletion) 


NW  IG(NW)  VTD  (NW)  J 

1  2  1 

2 

2  3  1 

2 

3  3  1 


2 


NI-1 

SIGD(NW,J,NI) 
NI=2  NI*3 

0.0 

0.0 

-  .105 

0.0 

0.0 

-1.32 

0.0 

-  .274 

-  .052 

0.0 

-1.560 

-  .660 

0.0 

-  .274 

-  .052 

0.0 

-1.560 

-  .660 
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controlled  by  one  of  the  control  programs  depending  on  the  WADOP  option 
chosen  (initialize,  add,  or  delete).  The  following  three  parts  each 
illustrate  the  operation  of  one  of  these  control  programs.  Each  control 
routine  utilizes  a  number  of  other  local  subroutines,  as  well  as  exter¬ 
nal  routines  (see  figure  56) . 

The  Add  Weapon  Control  routine  (Part  III)  will  be  used  as  a  vehicle  to 
illustrate  the  operation  of  the  program.  Once  the  operations  of  this 
routine  are  understood,  the  corresponding  operations  in  the  other  rou¬ 
tines  should  be  obvious. 

The  routine  first  checks  to  be  sure  that  the  number  of  weapon  additions 
and  weapon  deletion  operations,  IOP,  on  this  target  does  not  exceed  100. 
If  it  does,  it  is  assumed  that  STALL  and  WAD  are  caught  in  an  endless 
loop  probably  repeatedly  adding  and  deleting  the  same  weapon,  so  the 
processing  of  the  target  is  terminated.  In  principle,  such  looping 
should  not  occur.  However,  it  has  been  found  that  errors  in  reading 
file  data  for  one  target,  or  a  random  machine  malfunction,  or  inconsis¬ 
tencies  in  the  data  supplied  to  the  program  can  sometimes  result  in  such 
a  situation.  This  makes  it  possible  for  the  program  to  proceed  to  the 
next  target  rather  than  aborting  the  entire  run.  However,  when  this 
happens,  the  LOOP  flag  is  set  nonzero.  This  causes  the  print  in  state¬ 
ment  41  to  appear  during  the  initialization  of  every  succeeding  target, 
so  that  the  user  is  sure  to  notice  that  the  difficulty  occurred. 

However,  assuming  that  no  such  loop  has  occurred,  the  routine  adds  the 
Lagrange  multiplier  for  the  weapon  added  to  the  COST  of  the  allocation 
and  also  updates  SUMFREM,  the  sum  of  the  premiums  for  the  target.  (This 
last  variable  SUMPREM,  as  well  as  the  variables  TBENEFIT  and  TPMX  near 
statement  14,  are  computed  only  to  provide  a  consistency  check  on  the 
treatment  of  premiums .  These  variables  are  not  essential  to  the  opera¬ 
tion  of  the  program.)  Notice  that  these  variables  and  almost  all  other 
variables  used  by  WAD  such  as  VT,  VTP,  VTD,  PAYOFF,  PROFIT,  etc.  are 
computed  (even  for  multiple  targets)  as  if  the  program  were  dealing 
only  with  a  single  simple  target.  It  is  necessary  to  take  target  multi¬ 
plicity  into  account  only  when  dealing  with  variables  which  accumulate 
the  total  cost,  total  payoff,  or  total  consumption  of  specific  weapon 
groups  over  the  whole  target  system. 

The  PREMIUMS  used  by  the  allocator,  however,  depend  on  just  such  a  vari¬ 
able  —  namely  SURPWP(G),  the  available  surplus  (positive  or  negative) 
of  unused  weapons  in  each  group  G.  Consequently,  when  a  weapon  G  is 
added,  the  PREMIUM  for  that  weapon  group  must  be  recalculated. 

SURPWP(G)  -  SURPWP(G)  -1 
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The  program  proceeds  as  usual  to  update  SURPWP  for  the  weapon  now  being 
added  and  PREMIUMS (G)  is  called  as  usual. 

Both  the  add  and  delete  weapon  processes  exit  WAD  through  statement  12. 

At  this  statement,  WAD  calls:  1)  ADDSAL  to  update  the  salvoed  weapon 
stockpile;  2)  NEWSAL  to  determine  if  the  preferred  salvo  is  still  avail¬ 
able,  and  3)  WADOUT  to  calculate  the  decision  variables  for  STALL.  After 
these  calls,  WAD  updates  the  payoff,  cost,  profit  and  benefit  variables. 


m 


♦Local  subroutines 


Figure  56. 


Fart  III:  Add  Weapon  Control 
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18 


9 


Subtract  LAMEF(G)  From  Cost] 
For  Target,  And  Add  Premium) 
For  Deleting  (G)  To  Sum  Of 
Premiums  For  Target 


I 


Is  \  Yes 
PROGRESS. EQ.l? 


I 


No 


19_ 

For  Weapon  Group  Deleted 
Increase  SURPWP(G) 


920 


Call  PREMIUMS 
[Updates  PREMIUM(C)  To 
Reflect  Change  In 
SURPWP(G) 


I 


Call  CALSG* 

For  Each  Weapon  Group 
Calculates  Reduction  DSIG 
In  Correlation  Term 


20 


I 


Call  SUBS1G 
Updates  MU  And  SIG  In 
Payoff  Estimate,  Subtract 
DSIG  As  Required  In 
All  Estimates 


I 


Call  SUBIND* 
Updates  Indices  To 
Reflect  Deletion 
Of  Weapon  (G) 


22 


1 


Call: 

CALPAY" Actual  Payoff 
CALPOT  Weapon  Added  Payoffs 
CALDEL  Weapon  Deleted  Payoffsj 


♦Local  subroutine 


Figure  56.  Part  IV:  Delete  Weapon  Control 


292 


CH-2 


293 


3.11.9  Subroutine  WADOUT 


PURPOSE :  This  routine  summarizes  decision  alternatives  for 

STALL  by  combining  payoff  data  produced  by  WAD 
with  the  weapon  cost  data  (Lagrange  multipliers 
minus  premiums).  It  also  contributes  to  the  effi¬ 
ciency  of  WAD  by  making  inappropriate  weapons  in¬ 
active. 


ENTRY  POINTS: 
FORMAL  PARAMETERS: 
COMMON  BLOCKS: 

SUBROUTINES  CALLED: 
CALLED  BY: 


WADOUT 


None 

C30,  C33,  CONTRO,  DYNAMI,  PREMS,  SALVO,  WADFIN, 
WADOTX ,  WADWPN ,  WPFIX 

PRNTALL 

WAD 


Method : 

The  output  data  produced  by  WADOUT  consist  of  the  following  parameters 
for  STALL  that  are  recorded  in  /WADOTX/: 

PPMX  and  IPPMX  -  the  maximum  potential  increase  in  effective 
profit  for  any  single  weapon;  and  the  index  G  to  that  weapon 
group,  respectively. 

PVRMX  and  IPVRMX  -  the  maximum  effective  efficiency  of  any 
potential  weapon;  and  the  index  G  to  that  weapon  group,  re¬ 
spectively. 

DRfN  and  IDPMN  -  the  minimum  marginal  effective  profit  for  any 
weapon  now  assigned;  and  the  index  NW  to  that  weapon  in  the 
list  of  weapons  assigned,  respectively. 

It  also  produces  the  array  INACTIVE (G)  in  /WADWPN/  which  is  used  by  WAD 
to  determine  which  weapons  groups  need  not  be  processed. 

The  input  data  from  WAD  consist  of: 

VT  -  the  surviving  target  value  in  /C33/ 

VTD (N)  -  the  potential  weapon  deleted  surviving  target  value 
(also  in  /DYNAMI/  (equivalenced  to  RVAL) 

VTP(G)  -  the  potential  weapon-added  surviving  target  value  in 
/WADFIN/. 
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The  input  data  on  weapon  costs  consist  of  LAMEF(G),  PREMIUM(G),  and 
DPREMIUM(G)  . 

WADOUT  also  initializes  VTMAX  and  VTMIN  of  /WADOTX/ ,  and  MAXCOST  which 
reflect  the  MINKILL,  MAXKILL  specifications  for  each  target. 

The  quantities  ALPHA  and  VTEF  of  /WADOTX/,  are  essentially  local  vari¬ 
ables  for  WADOUT.  They  are  included  in  this  common  block  for  use  by 
PRNTNOW,  and  in  the  case  of  ALPHA,  to  allow  WAD  to  reinitialize  it  to 
1.0  for  each  new  target. 

The  flow  diagram,  figure  57,  is  in  three  parts.  In  part  I,  the  user- 
input  parameter  IMATCH  is  used  to  determine  the  method  of  computing 
MINKILL  and  MAXKILL.  If  IMATCH  is  0,  then  the  damage  calculations  used 
to  determine  residual  target  value  for  purposes  of  MINKILL  and  MAXKILL. 
use  time  dependence  of  target  value.  If  IMATCH  is  nonzero,  then  MINKILL 
and  MAXKILL  are  computed  relative  to  the  original  target  value.  In  addi¬ 
tion,  if  IMATCH  is  100,  then  the  routine  prints  the  WADOUT  variables,  VTO, 
VT,  VTZO  (original  target  value),  FLGMN ,  FLGMX,  SVTMIN,  SVTMAX,  VTMIN, 
VTMAX,  and  ALPHA. 

The  Do  loop  ending  at  statement  18  is  used  to  flag  all  groups  with  fixed 
weapons  on  this  target  as  active.  This  prevents  their  being  removed  dur¬ 
ing  processing  and  maintains  the  validity  of  the  damage  calculations. 

In  Part  II,  the  processing  begins  by  evaluating  the  actual  effective  sur¬ 
viving  target  value  VTEF.  It  then  scans  all  weapons  currently  assigned 
to  calculate  the  output  quantities  DMI  and  IDPMN.  If  any  weapon  now  on 
target  fails  to  destroy  a  fraction  of  the  original  value  greater  than 
MINDAMAG,  the  weapon  is  flagged  for  inmediate  removal  (statement  15). 

At  the  same  time,  the  groups  already  assigned  are  flagged  with  INACTIVE 
*  -100  to  eliminate  any  possibility  that  they  would  be  erroneously  set 
inactive.  (WADOUT  never  exits  with  INACTIVE  set  negative.  A  weapon 
group  flagged  with  a  -100  is  always  reset  to  INACTIVE  =  0  before  the  rou¬ 
tine  exits.)  (Do  14  loop) 

Basically,  the  processing  in  part  III  is  concerned  with  scanning  all  po¬ 
tential  weapons  to  calculate  the  output  quantities  PFMX,  IPPMX,  PVRMX, 
and  IPVRMX.  Any  weapon  which  would  fail  to  destroy  a  fraction  of  the 
original  value  greater  than  MINDAMAG  will  be  ignored  in  these  calcula¬ 
tions.  Thus,  it  could  never  be  allocated  to  the  target. 

In  addition,  if  a  weapon  is  a  salvoed  missile  but  no  salvo  is  available 
(i.e.,  MTSAL(G)  less  than  or  equal  to  zero),  or  if  the  group  has  no  weap¬ 
ons,  then  it  is  not  considered  as  a  potential  weapon. 

At  the  same  time,  however,  the  values  for  the  array  INACTIVE (G)  are  es¬ 
tablished.  The  INACTIVE  array  for  each  target  is  permanently  stored  on 
the  WPNTGT  files,  as  it  was  originally  computed  by  GETDTA,  with  only 
tvo  values  —  zero  for  weapons  in  range  of  the  target,  and  100  for 
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Figure  58.  Part  II:  (Part  2  of  2) 
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3.12.1  Subroutine  PRNTOD 

PURPOSE :  To  produce  optional  prints  for  overlay  DEFAL. 

(Options  27  and  28) 

ENTRY  POINTS:  PRNTOD 

FORMAL  PARAMETERS :  IOPT  -  Print  option  number 

COMMON  BLOCKS:  C30,  DEFCOM,  DEFRES 

SUBROUTINES  CALLED:  None 

CALLED  BY:  PRNTNOW 

Method: 

The  formal  parameter  IOPT  determines  whether  option  27  or  28  appears. 
The  result  of  the  option  appears  in  the  Users  Manual,  UM  9-77,  Volume 
III. 

Subroutine  PRNTOD  is  illustrated  in  figure  59. 
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Figure  63.  (Part  5  of  7) 
348 


Figure  66.  (Part  2  of  4) 


CH-2 


Figure  66.  (pert  3  of  4) 
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SECTION  5.  MODULE  ALOCOUT 


5-1  Purpose 

Module  ALOCOUT  is  responsible  for  selecting  optimum  DGZs  (desired  ground 
zeros),  also  called  weapon  aim  points,  for  weapons  allocated  to  target 
complexes.  ALOCOUT  also  resorts  weapon  assignments  at  the  group  level 
for  use  within  the  Sortie  Generation  subsystem.  If  executed  as  part  of 
mini-allocator  process  (DATAMAKE) ,  produces  a  STRIKE  format  tape  rather 
than  resorting  weapon  assignments.  (For  format  of  STRIKE  tape  see 
CSM  MM  9-77  Volume  IV). 

Module  ALOC  specifies  weapon  groups  assigned  to  targets  together  with 
associated  targeting  data.  ALOCOUT  extracts  data  from  these  records 
and  computes  any  aiming  offsets  required  by  the  plan.  For  simple  targets, 
no  calculations  are  performed.  In  the  case  of  complex  targets  which  can 
have  several  elements  at  slightly  different  coordinates,  ALOCOUT  employs 
subroutine  DGZ  (desired  ground  zero  selector)  to  select  optimum  aim  points 
within  the  target  complex. 

5.2  Input 


ALOCOUT  operates  after  module  ALOC  assigned  weapons  to  targets.  These 
records  (ASSIGN)  in  addition  to  the  supporting  data  base  structure  must 
be  defined  for  proper  execution. 

5.3  Output 

No  new  data  base  records  are  created  during  the  execution  of  ALOCOUT. 
However,  the  weapon  assignment  records  (called  ASSIGN)  are  modified  in 
two  ways.  First,  for  assignments  to  complex  targets  or  for  assignments 
to  cities  with  nonzero  radius,  offsets  as  determined  within  the  module 
are  included  within  the  ASSIGN  record.  Second,  the  assignment  records 
at  the  weapon  group  level  are  resorted  for  use  within  the  Sortie  Genera¬ 
tion  subsystem.  For  missile  groups,  the  sort  is  based  on  decreasing 
values  of  attribute  RVAL.  For  bomber  groups,  the  order  is  based  on 
penetration  corridor  index  and  within  the  corridor  sorted  based  on 
attribute  RVAL.  The  penetration  corridor  that  contains  the  largest 
number  of  strikes  appears  first  within  the  sort,  followed  by  the  pene¬ 
tration  corridor  of  the  next  largest  number  of  strikes  and  so  on. 

5.4  Concept  of  Operation 

ALOCOUT  (that  is,  subroutine  ENTMOD)  operates  with  three  overlay  links. 

The  first  overlay  reads  the  target  list  (TARNUM)  passes  controls  to  sub¬ 
routine  PROCCOMP  for  offset  calculations  when  applicable  and  finally 
supplies  optional  prints.  After  all  targets  have  been  processed  con¬ 
trols  passes  to  the  second  overlay  which  consists  entirely  of  subroutine 
SUMPRN  which  reorders  strikes  at  the  weapon  group  level  and,  if  requested, 
produces  prints  concerning  the  individual  assignments.  If  part  of  mini¬ 
allocator,  the  second  overlay  is  skipped  and  the  third  overlay  is  executed. 
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5.5  Identification  of  Subroutine  Functions 


5.5.1  Subroutine  PROCCOMP.  This  subroutine  controls  the  bulk  of  pro¬ 
cessing  for  offset  determination.  It  is  executed  by  subroutine  ENTMOD 
only  for  those  individual  targets  that  require  offset  calculations. 

After  offsets  have  been  determined  the  assignment  record  (ASSIGN)  is 
updated  to  include  the  values.  Then,  PROCCOMP  returns  to  subroutine 
ENTMOD  for  acquisition  of  the  next  target  and  the  associated  strikes. 

5.5.2  Subroutine  SUMPRN.  This  subroutine  constitutes  the  entire  second 
overlay  of  ALOCOUT.  Its  purpose  consists  of  resorting  the  weapon  strikes 
at  the  group  level  and  providing  optional  prints. 

5.5.3  Subroutine  MI  NX OUT.  This  subroutine  controls  the  third  overlay 
and  produces  a  STRIKE  tape. 

5 . 6  Common  Block  Definition 


Common  blocks  used  by  EVALALOC  are  outlined  in  table  11.  Common  blocks 
that  communicate  with  the  COP  are  given  in  appendix  A  of  Program  Main¬ 
tenance  Manual,  Volume  I. 


Table  11. 

ALOCOUT  Common  Blocks 
(Part  1  of  2) 

BLOCK 

VARIABLE  OR  ARRAY 

DESCRIPTION 

CITY 

ICITY 

Set  to  nonzero  for  targets  with  attri¬ 
bute  RADIUS  not  equal  to  zero 

Cl 

XO(J),  YO(J) 

Coordinates  of  target  element  J 

VI(J) 

Initial  target  element  values 

RADL(J) 

Lethal  radius  of  target  element  J 

VTOA(J, I) 

Value  of  target  element  J  immediately 
following  arrival  of  weapon  I 

S(J,I) 

Survival  probability  of  target  element 
J  relative  to  weapon  I 

VEFF(J.I) 

Effective  value  of  target  element  J 
relative  to  weapon  I 

X(I),Y(I) 

Offset  coordinates  of  weapon  I 

PDEL(I) 

Probability  of  delivery  of  weapon  I 

ERDEL(I) 

Error  in  delivery  of  weapon  I 

YDSCL(I) 

Scaled  yield  for  weapon  I 

VESC(I) 

Intermediate  computational  value  used 
in  subroutine  VAL  for  determination 
of  total  escaping  target  value 

NI 

Number  of  weapons  for  complex 

NJ 

Number  of  target  elements  for  complex 

GAMETIME 

KDAY 

Day  set  for  STRIKE  times  (set  to  1) 

KMON 

Month  set  for  STRIKE  times  (set  to  1) 

KYEAR 

Year  set  for  STRIKE  times  (set  to  1) 

HHR 

Hour  set  for  STRIKE  times  (set  to  10) 

GRPY 

GROUPY (250) 

Logical  switch  indicates  if  group  alio- 

cated 
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Table  11.  (Part  2  of  2) 


BLOCK  VARIABLE  OR  ARRAY 

IONPRT  IPINDAT 

PRINCE (9) 

ISKIPDGZ  ISKIPDGZ 


JAZ 

F  (400, 26) 

LOCFIN 

LOCFIN 

STRIKE 

TOA(I) 

IREFSTRK(I) 

N 

WPGT 

YDMIN 

IGRP 

DESCRIPTION 

User  supplied  print  frequency  for 
print  option  1 

Set  TRUE  if  user  requested  option 

Use  indicator  for  DGZ.  Normally  it- 
is  0.  Compress  resets  it  to  1  if 
more  than  50  calls  to  it  are  made  to 
reduce  the  number  of  target  elements 
for  a  complex  target;  DGZ  is  not  used 
again  for  the  target  in  this  case 

Holding  arrays  for  sort  purposes 

Starting  location  into  IRSET's  arrays 
for  adverb  FINDMIN  instructions 

Weapon  time  of  arrival  to  target 

Reference  code  of  weapon  strike 

Number  of  strikes 

Minimum  allowable  weapon  group  yield 
Group  number 
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Figure  69.  Subroutine  COMPRESS 
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5.7.2  Function  CUMINV 


PURPOSE: 


To  determine  the  value  X  such  that  Z  is  the 
probability  that  x  £  X. 


ENTRY  POINTS: 

FORMAL  PARAMETERS : 
COMMON  BLOCKS: 
SUBROUTINES  CALLED: 
CALLED  BY: 

Method: 


CUMINV 

Z  -  The  probability  that  x  s  x 

None 

None 

PROCCOMP 


Function  CUMINV  is  illustrated  in  figure  70.  By  definition. 


Z  =  P[x  <  X] 


2 


•  dt  for  0  <  Z  <  1 


CIMINV  uses  the  following  approximation  X  for  X: 


(  A!  ♦  A2  •  V  ♦  A3  •  V2 _ 

1  ♦  Bx  .  V  +  B2  •  V2  +  B3  .  V3 


where 


V  -  /  In  (W)  , 

Aj  *  2.S15S17 
A2  *  . 802S53 

A3  «  .010328 


Q  =  Z  or  1  -  Z  such  that  0  <  Q  < 


Bj  =  1.43278S 
B2  -  . 180269 
B3  =  . 00130S 


.5 


Initialize  XG  Array 
With  Value*  Of  Offset 
Aia  Point  Coordinates 
Of  Weapons 


1004 


Yes 


1003 


l 


10 

No/"*" 


One  Target 
Element  Or  Less 


1011 


Use  FINDMIN? 


D 
) 


i  Tes 


Call  PERTBLD 
To  Perturb 
Weapon  Laydown 


1012 


(loi^— * 


Call  FINDMIN 
To  Optimize  Weapon 
Laydown  In  Such  a  Way 
That  Surviving  Target 
Value  la  Minimized 


Reformat  Weapon  Aim 
Faint  Coordinate*  For 
Tranafer  To  PROCCOMF 


Call 

SEECALC 


(w?) 


Call  SEEINPUT  To 
Print  Data  Input 
From  PR0CC0MP 

~  n 

Call  SEECALC 
To  Print 
Laydown  Results 


Figure  72.  (Part  2  of  3) 
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5.7.5  Subroutine  FINDMIN 

PITRPOSE:  This  subroutine  uses  a  steepest  descent  method 

to  determine  a  local  minimum  of  a  function  of 
several  variables.  An  initial  estimate  of  the 
minimum  position  is  input,  together  with  various 
tolerances.  FINDMIN  uses  two  auxiliary  routines, 
F2BMIN  and  GRADF,  which  define  the  function  to  be 
minimized  and  its  gradient,  respectively.  DGZ 
uses  FINDMIN  to  find  the  DGZs  for  complex  targets. 

ENTRY  POINTS:  FINDMIN 

FORMAL  PARAMETERS :  XO  -  Initial  guess  at  aim  point  offsets 

N  -  Length  of  XO  vector 

IMAX  -  Maximum  number  of  iteractions  for  FINDMIN 

E1,E2  -  Tolerances  for  the  minimization 
X  -  Best  aim  point  offsets  as  determined  by 
FINDMIN 

FI  -  Minimum  value  found  for  escaped  target 
value 

IFLAG  -  Print  control  flag 

>0  :  FINDMIN  debug  prints  will  be  produced 

SUBROUTINES  CALLED:  F2BMIN,  GRADF,  SEECALC 


CALLED  B’S 


DGZ 


Method : 

Given  a  function  F(X^, 
an  initial  guess  (XO^, 


dF.  dF 

Xj),  the  gradients  G^  =  dX^  ,  and  G 2  *  ,  and 

XD^)  at  the  aim  point  offsets,  FINDMIN  finds  the 


local  escaped  target  value  F  and  its  associated  aim  point  offset  coor¬ 


dinates  (X^,  X2).  Each  iteration  consists  of  a  function,  F,  and  gradient, 
(G^,  G2) ,  evaluation  followed  by  determination  of  the  minimum  function 
value  along  the  line  associated  with  the  modified  steepest  descent  direc¬ 


tion.  F  is  redetermined  at  each  iteration  and  is  defined  in  such  a  way 


that  it  converges  after  two  iterations.  FINDMIN  uses  two  subroutines 
during  its  processing.  The  first,  F2BMIN,  defines  the  escaped  target 
value  function  F  in  terms  of  aim  point  offsets  X^  and  X2.  The  second, 
GRADF,  defines  the  gradient  components  G^  and  G2. 
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Subroutine  SEECALC  is  called  after  each  iteration  to  print  the  results 
of  the  optimization. 

Subroutine  F1NDMIN  is  illustrated  in  figure  74. 
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Figure  74.  (Part  2  of  4) 
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Figure  74.  (Pert  3  of  4) 
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Figure  74.  (Part  4  of  4) 
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5.7.10  Subroutine  PROCCOMP 

PURPOSE:  To  set  up  arrays  in  common  block  /Cl/  for  the  com- 

plex  target  so  that  Che  subroutine  DGZ  can  use  the 
arrays  during  the  selection  of  optimal  aim  point 
offsets  for  the  weapons  allocated  to  the  target; 
and  to  modify  target  weapon  assignments  records 
for  inclusion  of  the  computed  offsets. 

ENTRY  POINTS:  PROCCOKP 

FORMAL  PARAMETERS:  None 

COMMON  BLOCKS:  CITY,  Cl,  CIO,  C30,  ISKIPD,  STRIKE,  WPGT 

SUBROUTINES  CALLED:  COMPRESS,  CUMINV,  DGZ,  DIRECT,  ERGOT1,  ERGOT2, 

HEAD,  MODFY,  NEXTTT,  ORDER,  REORDER,  TIMEME, 

VALTAR 

CALLED  BY:  ENTMOD 

Method : 

When  ENTMOD  encounters  a  complex  target,  PROCCOMP  is  called  in  order  to 
assemble  data  in  a  form  that  can  be  efficiently  used  for  DGZ  selection. 
Each  target  component  of  the  complex  generates  a  standardized  "target 
element"  In  the  working  arrays  used  by  subroutine  DGZ  (common  /Cl/). 
Targets  with  more  than  one  hardness  component  generate  more  than  one 
such  target  element,  and  targets  with  a  specified  target  radius  will 
generate  several  elements  spread  over  the  area  of  the  target  to  repre¬ 
sent  a  value  over  the  area.  For  complexes,  individual  target  elements 
are  obtained  by  walking  the  data  base  chain  called  'CMPTGT'. 

If  the  number  of  target  elements  so  generated  exceeds  the  maximum  pro¬ 
gram  dimensions  (290),  subroutine  COMPRESS  is  called  to  recombine  target 
elements  near  each  other  having  nearly  the  same  lethal  radius.  In  any 
case,  for  efficiency  in  DGZ,  a  call  to  COMPRESS  is  made  just  before 
calling  DGZ.  On  return  from  DGZ,  PROCCOMP  modifies  weapon  assignment 
records  (ASSIGN)  for  definition  of  the  computed  offsets. 

Subroutine  PROCCOMP  is  illustrated  in  figure  79. 
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Figure  79.  (Fart  2  of  4) 
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5.7.11  Subroutine  SEECALC 

PURPOSE:  To  print  the  computation  values  relevant  to  the 

selection  of  aim  point  offsets  at  various  points 
within  the  DGZSEL  subarea  of  program  ALOCOUT. 

ENTRY  POINTS:  SEECALC 

FORMAL  PARAMETERS:  VESCIOT  :  Total  escaping  target  value 

XX  :  Vector  containing  the  aim  point  offset 

positions  for  the  weapons 

COMMON  BLOCKS:  Cl,  WAROUT 

CAT.T.CTI  BY.  DGZ ,  FINDMIN 

Method : 

When  called  by  DGZ  or  FINDMIN  subroutine  SEECALC  prints  the  title 
DGZSEL  COMPUTATION  VALUES  and  column  headings.  Then  for  each  weapon 
allocated  to  the  target,  SEECALC  prints  the  internal  weapon  number,  the 
aim  point  offsets,  and  the  survival  of  each  target  element  relative  to 
the  weapon.  At  the  end  of  the  print  for  each  target,  the  total  escaped 
target  value  is  printed. 

Subroutine  SEECALC  is  illustrated  in  figure  80. 
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Figure  80.  Subroutine  SEECALC 
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5.7.11.1  Subroutine  SEEINPUT 


PURPOSE: 

ENTRY  POINTS: 

FORMAL  PARAMETERS: 
COMMON  BLOCKS: 
SUBROUTINES  CALLED: 
CALLED  BY: 


To  print  target  data 
requested  by  ONPRINT 

SEEINPUT 

None 

Cl,  IONPRT 

None 

DGZ 


input  to  DGZSEL  when  so 
5. 


Method : 

If  ONPRINTS  5  was  specified  by  the  user,  SEEINPUT  prints  the  title  TARGET 
DATA  INPUT  TO  DGZSEL  DEBUG  PRINT  and  various  column  headers.  Subsequently, 
it  prints  data  from  iN  Cl. 

Subroutine  SEEINPUT  is  illustrated  in  figure  80.1. 
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Figure  80,1.  Subroutine  SEEINPUT 
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5.7.12  Subroutine  VAL 


PURPOSE: 


ENTRY  POINTS : 


VAL  determines  the  target  value  which  has 
escaped  for  a  given  weapon  configuration  and 
also  determines  the  effective  value,  F  ,  for 
each  target  element  as  seen  by  each  weapon. 

VAL 


FORMAL  PARAMETERS:  VESCTOT 


COMMON  BLOCKS:  Cl 


SUBROUTINES  CALLED:  None 

CALLED  BY:  DGZ,  F2BMIN 


Method : 

This  computation  uses  the  effective  values,  VEFF(J,I),  the  survival 
probabilities,  S(J, I),  and  the  time  dependent  target  values. 

Subroutine  VAL  is  illustrated  in  figure  81. 
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Figure  83.  Subroutine  WEPGET 
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5.8  Subroutine  SUMPRN 


PURPOSE: 


To  sort  weapon  group  assignment  chains  and  print 
out  assignment  summaries 


ENTRY  POINTS: 

FORMAL  PARAMETERS: 
COMMON  BLOCKS: 
SUBROUTINES  CALLED: 

CALLED  BY: 


SUMPRN 

None 

CIO,  C30,  GRPY,  IONPRT,  JAZ 

DIRECT,  DLETE ,  HEAD,  NEXTTT ,  ORDER,  REORDER, 
SORTIT,  TIMEME 

ENTMOD 


Method: 

This  subroutine  replaces  existing  ASSIGN  records  which  are  in  no  parti¬ 
cular  order  on  the  MYASGN  chain  with  the  same  set  of  ASSIGN  records  in 
the  following  sorts:  For  missile  groups,  assignments  are  sorted  on 
salvo  number  ascending,  and,  within  salvo,  on  the  value  of  the  RVAL 
attribute  descending.  For  bomber  groups,  assignments  are  sorted  on  cor¬ 
ridor  with  the  corridor  most  often  assigned  occurring  first,  and,  within 
corridor,  on  the  value  of  the  RVAL  attribute  descending. 

The  method  used  is  to  cycle  the  weapon  group  chain  and  perform  essen¬ 
tially  the  same  process  for  each  group.  First  a  record  is  read  from 
random  access  file  25.  This  record  contains  counts  of  the  weapon  groups 
assignments.  For  missile  groups  the  total  number  is  in  CORCNT(l)  and 
-1  in  C0RCNT(2) .  For  bombers  the  contents  of  CORCNT(I)  corresponds  to 
the  number  of  assignments  to  corridor  I.  At  this  point,  if  the  group  is 
a  missile  group,  the  print  header  is  produced.  If  the  group  is  a  bomber 
group,  the  assignments  are  totaled  and  the  proper  corridor  order  is 
stored  in  CORORD. 

Now  the  assignments  are  read  from  the  MYASGN  chain  and  stored  in  array 
F.  When  the  assignments  have  all  been  processed,  the  sort  keys  are 
generated.  Then  the  assignments  are  sorted  by  routines  ORDER  and  REORDER. 
The  old  assignments  are  deleted  from  the  MYASGN  chain.  Now  each  assign¬ 
ment  in  sort  order,  is  retrieved  from  array  F.  A  new  ASSIGN  record  is 
stored  in  the  MYASGN  chain  and  the  assignment  is  printed.  The  bomber 
header  is  produced  each  time  a  new  corridor  is  encountered. 

Subroutine  SUMPRN  is  illustrated  in  figure  84. 
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Figure  84.  (Part  3  of  12) 
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Figure  84.  (Part  8  of  12) 


Figure  84.  (Part  11  of  12) 


5.9  Subroutine  MINIOUT 


PURPOSE: 

Produce  miniallocator  STRIKE  tape 

ENTRY  POINTS: 

MINIOUT 

FORMAL  PARAMETERS: 

None 

COMMON  BLOCKS: 

CIO,  C15,  C30,  GAMETIME,  IONPRT 

SUBROUTINES  CALLED: 

ABORT,  CONVLL ,  DISTF,  FINDTIME, 
IGETHOB,  INFORM,  IPROB,  NEXTTT, 

HDFND,  HEAD, 
RETRV 

CALLED  BY: 

ENTMOD  (ALOCOUT) 

Method: 

This  subroutine  first  retrieves  the  necessary  weapon,  payload  and  warhead 
data.  Then  each  group  is  processed  in  order.  For  each  group,  every  asso¬ 
ciated  weapon  assignment  record  is  used  to  develop  and  reformat  the  neces¬ 
sary  data  to  produce  a  STRIKE  tape  record  (for  STRIKE  tape  format  consult 
CSM  MM  9-77  Volume  IV).  This  tape  is  thus  produced. 

Subroutine  MINIOUT  is  illustrated  in  figure  84.1. 


Main  subroutine  of  overlay  MINIO. 
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PINT-PKTX* ( ( PXLOW*RXLCM  ) + ( PXHIGH*RXHI GH  )  +  ( 1  - PXLOW- PXHIGH ) ]  *MI SDEF 

This  variable  is  Che  expected  number  of  objects  to  be  removed  by  the 
terminal  defense  interceptors. 

The  penetration  probability  for  any  warhead  is  defined  as 


If  this  probability  is  less  than  (1.0  -  FKTX)  it  is  reset  to  that  value. 
A. 2  Weapon/Target  Interaction 

The  quality  of  the  plans,  in  terms  of  realism  and  sophistication  ill 
be  a  direct  reflection  of  the  realism  Incorporated  in  the  payoff  junc¬ 
tion.  In  order  to  produce  plans  of  maximum  realism,  the  payoff  function 
should  reflect  all  the  major  factors  that  would  be  considered  by  an 
experienced  military  planner.  The  design  incorporates: 

1.  Time  of  arrival  of  weapons 

2.  Time  dependence  of  target  values,  which  can  reflect  a  planner's 
uncertainty  in  the  time  of  arrival  of  weapons  relative  to 
change  in  target  value 

3.  Weapon  range  limitations 

4.  Uncertainty  in  target  vulnerability 

5.  Correlations  in  the  effectiveness  of  weapons  of  similar  nature 
reflecting  such  factors  as  reliability,  DBL  probability,  and 


defense 

effectiveness. 

To  evaluate  the  capability  of  any  weapon  group  against  any  target,  pro¬ 
gram  ALOC  requires  six  basic  numbers.  These  are: 

SBL(G) 

The  probability  assumed  that  weapons  in  group  G  are 
not  destroyed  before  launch 

CC(KR) 

The  assumed  casnand  and  control  reliability  associ¬ 
ated  with  the  region  for  group  G 

REL(K) 

The  assumed  reliability  for  weapon  type  K  used  by 
group  G 

PEN(G) 

The  estimated  penetration  probability  for  weapons 
from  group  G  to  the  target 

STK(G, JH) 

The  estimated  kill  probability  of  warheads  in  group  G 
if  delivered  against  the  JH  hardness  component  of  the 
target 
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TVALTOA(G) 


The  estimated  target  value  at  the  time  of  weapon 
arrival  for  weapon  from  group  G  (this  factor  is  com¬ 
puted  from  the  time  of  arrival  for  a  weapon  from 
group  G,  TOA[G] ) . 

These  numbers  reflect  the  planning  factors  the  user  has  specified  for 
the  plan  generation  and  do  not  necessarily  reflect  the  values  that  the 
user  specifies  for  external  simulation.  The  number  is  noted  as  "assumed 
where  it  is  a  direct  user  input  supplied  in  the  data  base.  It  is  de¬ 
scribed  as  "estimated"  where  it  is  a  derived  quantity,  based  on  other 
input  data. 

Actually,  the  numbers  reflect  only  two  types  of  information  --  the  time 
of  arrival  information,  and  the  kill  probability  data.  The  single  shot 
kill  probability  is  simply  a  product  of  the  first  five  items.  The  break 
down  of  the  single  shot  kill  probability  into  these  five  separate  fact¬ 
ors,  however,  is  required  in  order  to  estimate  correlations  in  delivery 
probability  between  several  warheads  delivered  to  the  same  target. 

Most  of  the  processing  of  weapon/ target  interactions  deals  with  the  six 
quantities  given  above.  These  quantities  are  then  used  in  the  calcula¬ 
tion  of  weapon  payoff. 

The  basic  payoff  calculation  is  modified  by  the  inclusion  of  weapon  cor¬ 
relation  considerations.  For  each  single  weapon,  four  factors  are  cal¬ 
culated:  the  single  shot  kill  probability  and  three  auxiliary  quanti¬ 
ties  required  by  the  correlation  model  (see  Weapon  Correlations,  this 
appendix) . 

If  we  define  the  overall  single  shot  kill  probability  on  one  hardness 
component  J  as:  SSK  *  REL  *  CC  *  SBL  *  PEN  *  STK 
then  MUP(G.J)  -  -  L0GF(1.0  -  SSK) 

and  SSIG(G,J)  -  MUP(G,J)/-LOGF(SSK) 

If  the  option  to  use  the  square  root  damage  law  is  selected,  MUP  is 
defined  in  a  different  manner.  It  is  defined  so  that: 

(1.0  -  SSK)  -  (l  +yrMUP(G,J)j  *exp  VMUP(G,J)j 

The  use  of  the  square  root  damage  function  is  further  explained  in  a 
later  section  (see  Multiple  Weapon  Attacks  —  Square  Root  Law,  this 
appendix) . 

MUP  is  in  effect  a  measure  of  the  effectiveness  of  the  weapon  against 
the  specified  hardness  component.  If  all  weapons  were  independent,  the 
survival  probability  for  the  component  with  respect  to  multiple  weapons 
IG  would  be  simply: 

EXPF  -(E  MUP(IG, J)) 

(This  is  called  the  exponential  damage  law.) 
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where  MAXSAL  Is  the  maximum  salvo  number  in  the  group  and  .1  is  a  sensi¬ 
tivity  parameter  to  slow  the  rate  of  change  of  P. 

Closing  Factors  —  Premiums:  The  Lagrange  multiplier  for  each  weapon  is 
modified  by  a  premium.  This  factor  is  used  to  force  closure  of  weapon 
allocations  to  the  available  stockpile.  It  acts  as  a  bonus  for  using 
under-allocated  weapons  and  a  penalty  for  using  overallocated  weapons. 
The  parameters  which  are  used  to  calculate  the  premiums  are: 

SURPWP(G)  An  estimate  of  the  number  of  surplus  (or  unallocated 
weapons)  in  the  group.  This  number  is  based  on  esti¬ 
mated  allocation  rates  in  the  early  phase  and  the 
actual  allocation  later. 

NWPNS(G)  The  actual  number  of  weapons  in  group  G. 

LAMEF(G)  The  Lagrange  multiplier  for  the  group. 

The  premium  depends  also  on  three  control  parameters:  PROGRESS,  PRM, 
and  CLOSE.* 

The  effect  of  PROGRESS  (described  earlier)  is  as  follows; 

1.  If  PROGRESS  is  greater  than  1.0,  this  indicates  that  a  verifi¬ 
cation  allocation  is  desired  to  obtain  a  theoretical  upper 
bound  on  the  payoff  without  regard  to  meeting  the  actual  stock¬ 
pile  constraints.  For  this  purpose,  the  premiums  are  simply 
set  to  zero. 

2.  If  PROGRESS  is  less  than  1.0,  a  small  premium  is  computed  which 
is  intended  only  to  avoid  large  deviations  from  the  desired 
allocation  rate  of  small  errors  in  the  Lagrange  multipliers. 

(Otherwise,  a  trivial  change  in  the  multipliers  for  two  com¬ 
peting  weapons  could  result  in  a  complete  change  from  always 
allocating  one  to  always  allocating  the  other.) 

3.  If  PROGRESS  is  equal  to  1.0,  this  is  a  signal  that  the  closing 
phase  has  been  reached  and  the  object  is  to  close  in  on  an 
exact  allocation  of  the  available  weapons.  In  this  case,  a 
larger  step  function  premium  is  computed,  and  the  size  of  the 
step  function  is  gradually  increased  until  final  closure  occurs. 


PROGRESS  is  set  Internally  by  the  module  as  described  in  Section  3. 
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During  the  early  allocation  phase,  superimposed  on  the  actual  payoff  is 
a  small  negative  quantity  (called  a  premium)  that  is  proportional  to  the 
value  of  each  weapon  group  and  quadratic  in  the  size  of  the  error  in 
allocation.  In  effect,  the  actual  payoff,  H(X),  for  any  allocation,  X 
is  adjusted  to  H(X)ADj: 

12  j 

NWPNS  (G)*LAMEF  (G)*  (S^s(G)) 

) 

This  quadratic  addition  to  the  payoff  function  has  the  effect  of  intro¬ 
ducing  a  preference  for  allocations  where  the  absolute  value  of  SURPWP 
is  small. 

The  addition  or  deletion  of  a  weapon  from  group  G  will  give  rise  to  a 
difference  in  SURWP  equal  to  the  current  target  multiplicity.  Thus, 
the  change  in  this  quantity  (per  unit  multiplicity)  with  the  addition 
of  a  weapon  G  is: 


PREMIUM(G) ”PRM*LAMEF (G) * 


SURPWP (G)  -  .5 


NWPNS (G) 


and  the  change  with  deletion  of  a  weapon  is: 


DPREMIUM ( G) -PKM*LAMEF ( G) * 


NWPNS (G) 


The  value  of  PKM  is  a  user- input  parameter.  The  value  should  be  less 
than  1.0.  Otherwise,  in  cases  when  no  weapons  from  some  group  have  been 
used,  the  premium  for  allocation  of  a  weapon  could  exceed  the  cost  of 
the  weapon  LAMEF(G)  and  weapons  could  be  allocated  even  if  the  payoff 
were  zero  or  even  negative.  Experience  has  shown  that  values  between 
.5  and  .9  work  very  well. 

When  PROGRESS  reaches  1.0,  PKM  is  set  to  .9  by  the  program  to  accelerate 
convergence.  In  addition,  a  small  step  function  is  added. 

The  following  sketch  illustrates  the  value  of  these  step  function  pre¬ 
miums  as  a  function  of  their  SURPWP: 
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EXPASM  ”  fraction  of  weapons  in  a  group  which  are  ASMs 

EXPBMB  -  1  -EXPASM  •  fraction  of  weapons  which  are  bombs 

DEA  -  expected  destroyed  value  of  target  if  ASM  used  (calculated 
in  program  ALOC) 

DEB  *  expected  destroyed  value  of  target  if  bomb  used  (calcu¬ 
lated  in  program  ALOC) 

AVDE  ■  average  (by  group)  of  quantity  ABSF  (DEA  -DEB) 

FASM  *  current  fraction  of  weapons  allocated  which  are  ASMs  (cal¬ 
culated  in  program  ALOC) 

FBOMB  »  1  -FASM  *  current  fraction  of  weapons  allocated  which  are 
bombs 

CONPAY  “  internal  program  variable  between  0.  and  1.  ; 

Except  for  CONPAY,  all  of  the  above  variables  are  defined  for  each  group 
composed  of  bombers.  For  each  bomber  group  on  each  target  the  alloca¬ 
tion  process  selects  the  type  of  warhead  (ASM  or  bomb)  which  is  to  be 
used  on  the  target.  Uhen  the  value  of  PROGRESS  is  zero  or  two,  the  pre¬ 
ferred  weapon  is  the  weapon  with  the  higher  DE  (l.e.,  ASM  will  be  selec¬ 
ted  if  DEA  is  greater  than  DEB  and  vice  versa).  For  values  of  PROGRESS 
of  .4,  .5,  .75,  and  1.0,  the  selection  process  will  consider  the  alloca¬ 
tion  franctions  of  the  ASMs  and  bombs,  as  described  in  the  following 
paragraphs . 

If  ASMs  are  underallocated,  ASMs  are  selected  as  the  preferred  weapon 
unless  DEB  is  greater  than  DEA  and  the  quantity  (EXPASM  -FASM) /EXPASM 
is  less  than  or  equal  to  the  quantity  CONPAY* (DEB-DEA) /AVDE.  If  both 
of  these  two  conditions  are  met,  then  the  preferred  weapon  is  the  bomb. 

Note  that  the  quantity  (EXPASM  -FASM) /EXPASM  provides  a  measure  of  the 
size  of  the  allocation  imbalance.  If  ASMs  are  only  slightly  underallo¬ 
cated,  this  quantity  will  be  very  small  (near  zero).  If  there  is  a 
great  difference  between  the  actual  and  allocated  fractions,  this  quan¬ 
tity  will  approach  the  value  one.  The  quantity  (DEB-DEA) /AVDE  is  a  meas¬ 
ure  of  the  magnitude  of  the  damage  difference  relative  to  the  average 
damage  difference.  The  quantity  ranges  from  a  low  of  zero  to  high  posi¬ 
tive  values.  A  value  of  one  for  the  quantity  represents  an  average 
damage  difference.  The  variable  CONPAY  is  used  to  reflect  the  impor¬ 
tance  of  the  allocation  difference  relative  to  the  damage  difference. 
Thus,  the  conditions  of  the  preceding  paragraph  select  the  bomb  as  the 
preferred  weapon  if  the  allocation  difference  is  less  than  the  modified 
damage  difference.  Thus,  if  the  allocation  is  nearly  correct,  the  more 
damaging  weapon  is  likely  to  be  chosen  as  preferred.  If  the  allocation 
is  far  from  correct,  the  underallocated  weapon  will  be  selected  on  all 
targets  except  those  where  the  damage  difference  is  quite  large.  The 
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same  rationale  holds  for  the  case  of  underallocated  bombs  as  described 
in  the  next  paragraph. 

If  bombs  are  underallocated,  bombs  are  selected  as  the  preferred  weapon 
unless  DEA  is  greater  than  DEB  and  the  quantity  (EXPBMB  -FBOMB) /EXPBMB 
is  less  than  or  equal  to  the  quantity  CONPAY*(DEA  -DEB)/AVDE.  If  both 
these  conditions  are  met,  then  the  preferred  weapon  is  the  ASM. 

The  value  of  variable  CONPAY  lies  between  zero  and  one.  Lower  values  of 
CONPAY  tend  to  increase  the  importance  of  the  allocation  difference. 

High  values  of  CONPAY  increase  the  importance  of  the  damage  difference. 

In  order  to  provide  adequate  closing  force,  the  value  of  CONPAY  decreases 
as  the  value  of  PROGRESS  increases.  Additionally,  when  PROGRESS  equals 
one,  the  value  of  CONPAY  continues  to  decrease. 

The  selection  of  ASM  or  bomb  on  a  particular  target  allows  the  alloca¬ 
tion  process  to  assess  correctly  the  expected  damage  effects.  Bombs 
and  ASMs  usually  differ  greatly  in  yield,  penetration  probability,  CEP, 
and  delivery  probability.  By  differentiating  between  these  weapon  types 
at  allocation  time,  the  allocation  program  selects  the  best  weapon  to 
be  used  when  the  bomber  sorties  are  generated.  The  balance  between 
allocation  and  damage  differences  provides  for  maximization  of  damage 
while  continuing  consideration  of  actual  weapon  stockpiles. 

A. 5  Derivation  of  Lagrange  Multiplier  Adjustment 

Define  the  following  variables: 

*  sum  of  the  target  weights  from  last  Lagrange  mul¬ 
tiplier  update 

-  number  of  weapons  sharing  attribute  J 

■  number  of  targets 

-  user-input  parameters  which  control  rate  of  multi¬ 
plier  adjustment 

■  Lagrange  multiplier  for  group  G 

•  Lagrange  multiplier  for  attribute  J;  J  -  ALL' 

CLASS,  TYPE,  etc. 

-  local  Internal  control  variable  which  governs 
size  of  premiums  (closing  factors) 

■  number  of  weapons  in  group  G 


CURSUM(J) 

NOW  PS  (J) 

NTGTS 

SNSTVTY 

FSNSTVTY 

LAMEF(G) 

LA(J) 

PPM 

NWPNS(G) 
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If  we  now  associate  the  variable  PARTIAL  with  n  this  gives  rise  to  the 
following  procedure  for  updating  LA: 


T  a  i  T\  ,TA  (  n  *  M  n  4.  CORFAC*ALERREST(J.  INTPRD )/(-  PARTIAL), 
la0w;  li,u  +  ALERREST(J,  INTPRD)  +  (NCWPS(J) /NTGTS )J 


ALERREST(J)  is  computed  as 


ALERREST(J) 


RUNSUM(J)  NOWPS (  J) 

WTSUM(J)  "  NTGTS 


The  formula  for  LAi (J)  is  well-behaved  if  ALERREST  is  large  and  positive 
but  if  it  is  negative  and  as  large  as  the  expected  rate  (NOWPS (J) /NTGTS) 
(i.e.,  if  the  actual  allocation  rate  is  zero),  then  the  denominator  goes 
to  zero.  In  this  case  an  infinite  correction  would  be  indicated.  To 
avoid  this,  the  expected  rate  in  the  denominator  is  multiplied  by  2 
giving: 

-  ,  ...  ..  (nie  n  0  CORFAC*ALERREST ( J ,  INTPRD) / (-PARTIAL)  . 

1W'  ”  +  ALERREST(J,  INTPRD)  +  2* (NOWPS (J) /NTGTS )J 

This  is  the  function  used.  The  new  Lambda's  LA^(J),  are  recomputed  for 
attribute  J  (e.g..  Jail,  Jclass)  and  for  every  MULSTEP  targets  as  pre¬ 
viously  outlined. 

In  the  present  version  of  the  program  the  value  of  PARTIAL (J)  has  been 
set  equal  to  1.0  for  all  the  local  multipliers  LA(J).  This  choice  is 
based  on  the  effect  of  the  return  on  the  sensitivity  of  the  allocation 
rate  to  the  value  of  LAMEF  or  A.  When  the  multipliers  are  almost  cor¬ 
rect,  it  is  usually  the  case  that  most  weapon  groups  are  in  close  com¬ 
petition  with  many  other  groups  with  very  similar  properties.  Then  a 
small  change  in  the  multiplier  LAMEF  will  produce  a  very  large  change  in 
the  allocation  rates,  as  the  weapon  group  in  question  almost  totally 
replaces,  or  is  replaced  by,  its  competitors. 


However,  such  a  large  error  in  the  allocation  rate  will  not  actually 
occur  because  as  the  error  builds  up  the  estimated  value  of  the  payoff 
will  be  automatically  changed  by  the  premium.  Thus  for  constant  values 
of  LAMEF,  when  an  equilibrium  allocation  rate  is  reached,  it  must  be 
approximately  true  that  the  error  in  LAMEF  is  compensated  by  the  pre¬ 
mium.  This  is,  if  \q  is  the  correct  value  for  LAMEF  then: 

LAMEF  -PREMIUM  ~  AQ 


Since: 


PREMIUM 


PFM*LAMEF* 


SURPWP  -.5 
NWPNS 


we  can  define  a  relation  between  LAMEF  and  (SURPWP/ NWPNS) 
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LAMEF*(1  -PRM* 


;,URPWP  -,5v  ~  x 
NWPNS  ;  “  0 


Since  this  relationship  is  the  same  for  all  groups  it  is  reasonable  sim¬ 
ply  to  use  the  same  value  1.0  of  partial  derivative  for  all  local  multi¬ 
pliers. 


The  values  of  LAMEF(G),  where  G  is  the  group  index,  are  recomputed  using 
the  new  values  of  the  local  multipliers  (LA(J)  accordingly, 


LAMEF(G)  *  LA(Jall)*LA(Jclass)*LA(Jreg)*LA(Jalert)*LA(Jgroup) 

At  the  same  time  it  is  necessary  to  reevaluate  the  summation  of  the 
value  of  all  the  weapons  VALWPNS  ■  £  LAMEF(G)*NWPNS(G)  and  the  sumnation 
of  the  value  of  the  error  in  weapons  allocated 

VALERR  -  2  LAMEF(G)*ABSF (SURPWP(G) ) 

using  the  updated  values  of  LAMEF. 

Target  Weight  Change  Rate  &  Integration  Period:  The  above  explained 
how  multipliers  are  recomputed  by  monitoring  allocation  rates.  The 
remaining  discussion  addresses  how  the  target  weight  change  rate  and 
integration  period  is  computed. 

The  average  number  of  targets  over  which  allocation  rates  are  averaged 
(the  integration  period)  is  determined  by  the  rate  at  which  the  target 
weights  are  increased. 

In  estimating  the  rate  with  which  to  correct  multipliers,  it  was  com¬ 
puted  on  a  statistical  basis  that  even  if  the  allocation  rates  were 
correct  an  estimated  error  of  size  ALERREST  would  be  expected  if  the 
allocation  rates  were  monitored  only  over  a  small  sample  of  M  targets 
where : 

M  -  (NOWPS (J) /NTGTS) / (ALERREST(J))2 

Thus  if  separate  integration  periods  could  be  used  for  each  local  multi¬ 
plier,  M  as  defined  above  might  provide  a  reasonable  basis  for  deter¬ 
mining  the  period.  However,  in  fact,  the  same  periods  must  be  used  for 
all  local  multipliers  LA(J).  Currently  three  periods  are  maintained 
(INTPRO*!,  2,  3).  Consequently  the  value  of  the  integration  period  used 
must  be  based  on  an  estimate  of  overall  error  rate.  The  corresponding 
relation  is: 

M  -  (2  NOWPS (J) /NTGTS )/2  (ALERREST(J))2 
G  G 


where  the  simulations  are  taken  over  all  weapon  groups.  The  quantity 
ZNOWPS(J),  is  identical  with  NOWPS (1)  (Note:  LA(J)  for  J  -  I  is  used 
G 
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From  the  previous  equations 


fi  ■[&  sk  v<Ti)  “>d  Ei  -[j^ 


VCT.^)  -  V(T 


(SS) 1 


is  understood  =  1.  Also  V(TN  + 


(For  i  *  1,  the  product 
The  total  escaping  value  associated  with  target  j  is 
N-l  N+l 


0.) 


E^  »  E 

i«l  3  i«l 


/ 

’  i-1  ’ 

r  i 

11  Ski 

Lk-1 

|yw  -  vvj 

The  value  on  target  j  which  escapes  after  arrival  of  weapon  i  is  given 
by 

N+l 

P.?.i 


The  effective  value  of  target  j  associated  with  weapon  i  defined  by 


This  value  is  introduced  for  computational  efficiency  and  may  be  thought 
of  as  the  total  value  available  for  weapon  i,  the  effect  of  all  other 
weapons  having  been  taken  into  account. 

The  marginal  value  picked  up  on  target  j  due  to  weapon  i  is  given  by 


V1  -  V 


where  Sjj  is  a  function  of,  among  other  things,  the  position  of  weapon  i. 
For  a  fixed  weapon  configuration,  weapon  i  can  be  moved  from  (x,y)  to 
(x',y')  and  the  marginal  escaped  value  is  given  by: 

NT 

£  '  s'ij> 
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To  establish  an  initial  weapon  configuration,  a  lay-iiwn  is  performed  as 
follows.  Initially,  set  *  1  for  all  i,  j.  Denote  by  the  sur¬ 
vival  probability  of  the  k*h  target,  relative  to  the  ith  weapon,  when 
this  weapon  is  placed  on  the  target.  Now  the  ith  weapon  is  placed 
on  that  target,  j,  which  yields  a  maximum  value  for  the  expression 

NT 

£  Fik<sik  -  sik> 

The  are  now  set  to  equal  to  s4^  (k  ■  1,2 . NT)  the  (all  i,  k) 

are  redetermined,  i  is  increased  by  one,  and  the  process  repeated  until 
all  weapons  have  been  allocated. 

This  weapon  configuration  can  now  be  input  as  the  initial  position  to  a 
"hill  climber"  routine,  based  on  a  steepest  descent  algorithm,  which 
attempts  to  optimize  further  by  replacing  the  discrete  set  of  possible 
weapon  positions  with  the  two-dimensional  continuum.  The  function  to  be 
minimized  is: 


NT 

£ 


N+l 

£ 


j-1  i=l 


Eij 


Processing  by  the  optimizer  will  be  terminated  either  when  the  optimum 
has  been  achieved  or  when  a  specified  number  of  iterations  have  been 
completed.  In  either  case,  to  insure  that  the  local  optimum  obtained 
cannot  be  further  improved,  the  value  of  removing,  in  sequence,  each  of 
the  weapons  from  its  final  location  and  placing  it  on  one  of  the  target 
points  is  explored.  If  the  results  obtained  by  this  method  are  better 
than  those  achieved  with  the  previous  configuration,  this  new  assignment 
will  be  used  as  an  initial  one  for  a  second  utilization  of  subroutine 
FINDMIN.  If  not,  the  results  of  the  first  use  of  subroutine  FINDMIN 
will  be  kept. 
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